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11.  ABSTRACT 


The  impulse  response  of  a  kinetically  defined,  non-uniform  electron  plasma  slab 
bounded  by  sheaths  is  obtained.  The  one-dimensional  Vlasov-Poisson  system  is  consid¬ 
ered.  Electron-free  boundaries  are  specified  for  all  time.  Possible  equilibrium 
states  are  studied;  under  relatively  loose  mathematical  conditions  the  electron  dis¬ 
tribution  is  separable  and  the  only  possible  velocity  distribution  is  Gaussian.  Ad¬ 
missible  initial  conditions  are  investigated.  The  equilibrium  plasma  is  stimulated 
by  an  impulse  electric  field  which  imparts  upon  all  electrons  an  equal  additional  ve¬ 
locity.  A  numerical  scheme  is  developed  which  permits  to  obtain  reliable  solutions 
for  certain  time  periods  after  which  numerical  instabilities  begin  to  appear.  Througj>| 
out  this  study,  considering  plasma  slabs  with  an  equilibrium  temperature  around 
30000OK  and  a  plasma  frequency  of  approximately  1  GHe  at  the  center  of  the  slab,  re¬ 
liable  solutions  were  obtained  for  the  spans  ranging  from  2  plasma  periods  for  strong 
stimuli  to  20  plasma  periods  in  linear  operation. 

It.  is  found  that  for  excitations  in  the  linear  region  the  Fourier  Spectrum  of 
the  impulse  response  contains  three  strong  resonances  at  and  below  the  background 
plasma  frequency  at  the  center  of  the  slab  as  well  as  a  series  of  smaller  resonances 
above  this  frequency.  The  former  are  Tonks-Dattner  resonances  and  it.  is  shown  that 
each  one  is  associated  with  a  well-definable  region  within  the  slab.  The  highest  of 
the  three  high-amplitude  resonances  is  a  cold  plasma  oscillation  while  the  other  two 
are  shown  to  be  thermomodes  in  the  sense  of  Baldwin  and  Ignat.  In  the  non-linear 
responses  corresponding  to  strong  stimuli,  early  bunching  and  wave-particle  inter¬ 
action  are  observed. 
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ABE  "RACT 

The  impulse  response  of  a  kineiiically  defined,  non-uniform  elec¬ 
tron  plasma  slab  bounded  by  sheaths  is  obtained.  The  one-dimensional 
Vlasov-Foisson  system  is  considered.  Electron-free  boundaries  are  specified 
for  all  time.  Possible  equilibrium  states  are  studied;  under  relatively 
loose  mathematical  conditions  the  election  distribution  is  separable  and 

th;  only  possible  velocity  di stribution  i s  Gaussian.  Admissible  initial _ 

conditions  are  investigated.  The  equilibrium  plasma  is  stimulated  by  an 
impulse  electric  field  which  imparts  upon  all  electrons  an  equal  additional 
velocity.  A  numerical  scheme  is  developed  which  permits  to  obtain  reliable 
solutions  for  certain  time  periods  after  which  numerical  instabilities  begin 
to  appear.  Throughout  this  study,  considering  plasma  slabs  with  an  equilib¬ 
rium  temperature  around  300C0°K  and  a  plasma  frequency  of  approximately  1  GHs 
at  the  center  of  the  slab,  reliable  solutions  were  obtained  for  time  spans 
ranging  from  2  plasma  periods  for  strong  stimuli  to  20  plasma  periods  in 
linear  operation. 

It  is  found  that  for  excitations  in  the  linear  region  the  Fourier 
spectrum  of  the  impulse  response  contains  three  strong  resonances  at  and 
below  the  background  plasma  frequency  at  the  center  of  the  slab  as  well  as 
a  series  of  smaller  resonances  above  this  frequency.  The  former  are  Tonka- 
Dattner  resonances  and  it  is  shown  that  each  one  is  associated  with  a  well- 
definable  region  within  the  Blab.  The  highest  of  the  three  high-emplitude 
resonances  is  a  cold  plasma  oscillation  while  the  other  two  are  shown  to  be 
thermomodes  in  the  sense  of  Bladvin  ana  Ignat.  In  the  non-linear  responses 
corresponding  to  strong  stimuli,  early  bunching  and  wave-particle  interaction 


are  observed. 
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LIST  OP  SYMBOLS 


NOTE: 

1.  All  waits  are  expressed  In  the  MKS  International  system 
unless  noted. 

2.  Notation  in  the  appendices  is  independently  defined,  therein. 

3.  Indexing  superscripts  cad  subscripts  are  not  listed. 

a  Thar max  velocity  of  electrons  in  the  plasma  slab  in  equilibrium 

d  ifalf-width  of  slab 

D  S/3 t  (Operator) 

e  Electron  charge,  absolute  value 

E  Electric  field 

E(x)  Ditto,  in  the  equilibrium  state 

E(x,t)  Ditto,  time  dependent 

{E  }  Set  formed  by  E(x)  and  its  first  n  derivatives 

n 

**x 

f  In  Parts  IV  and  V:  Frequency  referenced  to 

i  In  Parts  I-III:  Electron  distribution  function 

f(x,v)  Ditto,  in  the  equilibrium  state 

f(x,v,t)  Ditto,  time-dependent 

f  Frequency,  referenced  to  T  * 

n  p 

fp  Plasma  frequency 

f  Ditto,  at  center  of  plasma  slab  in  equilibrium 

po 

h  Discrete  time  interval 

h(x)  Defining  function  of  !p(x),  in  order  that  tp(x)  can  be  used 

as  initial  condition  (See  Eq.  35) 

J  Total  number  of  time  intervals  covered  by  a  numerical  experiment 


viil 


m  Electros  mass 


M  Half-number  of  Intervals  covering  the  width  of  the  slob 

n,n(x)  Ion  density 

nfi  Electron  density  (See  also  <p(x)) 

n£(x,fc)  Ditto,  time-dependent 

-  H  Helf»nuabar  of  Intervals  in  (truncated)  velocity-space 

Total  number  of  electrons  per  unit  area  of  slab,  normalized 
to  unity  at  t“0 

N^C  H^-1.004  (See  text,  Section  Ill.b) 

Nm  Total  number  of  electrons  per  unit  area  of  slab 

* 

Nq  Electron  density  at  center  of  plasma  slab  in  equilibrium 

Q  ,  Coefficients  of  powers  of  w  (See  Eo.  28) 

nj 

R  n-th  partial  time  derivative  of  E(x,t)  at  t-0, 

n 

a  Order  of  infinity  of  the  voltage  on  the  boundaries  (walls) 


In  Darts  IV  and  Vs  Time  referenced  to  T 


In  Parts  I-III: 


T 

P 

u 


v 


V 

V 


V(x) 

V(t) 


Time  referenced  to  the  plasma  period 

Plasma  period  at  center  of  slab  in  equilibrium 

Fixed  value  of  the  velocity  imparted  to  all  electrons  by 
the  impulse  electric  field  at  t*0 

Electron  velocity 

Absolute  value  of  the  limits  of  truncated  velocity  space 

In  Ports  IV  and  V:  Identical  to  V  (See  ibidem.) 

n 

In  Parts  I-III:  Voltage 

Ditto,  in  equilibrium  state  V(o)»0 
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V  Voltage  normalized  to  the  parameters  of  the  standard 

u  experiment  (See  Eq.  62) 

w  v+u 

x  distance  from  center  of  slab 

x  .  First  minimum  of  ion  density  nearest  to  boundary  (wall)  x^d 
win 

v  Distance  from  wall 

Free  space  permittivity 

JyMx)  Debye  lanpth  in  equilibrium  state 

Ditto,  at  center  of  slab 
Do 

<p,(p(x)  Electron  density  in  the  plasma  slab  in  equilibrium  (See  also 
n  »  n  (x,t)) 

6  C 

{tpn>  Set  formed  by  <p(x)  and  its  first  n  derivatives 

^i(v)  Velocity  distribution  in  the  plasma  in  equilibrium,  normalised  to 

/^“’Kvjdv  «  1 
-00 

Plasma  angular  velocity 

u  Ditto  at  center  of  plasma  slab  in  equilibrium 

po 


I.  -  IMTRQDUCTIrtN  ANH  HISTORICAL  BACKGItOtJED 

The  study  of  plasmas  is  of  great  interest  because  it  covers  a 

wide  spectrum  a £  phenomena  which  ore  applicab3.e  to  many  fields  from 

fundamental  research  to  industrial  applications.  To  mention  a  few  areas 

of  investigation!  the  studies  of  plasmas  ranges  from  the  fcicrocoemie 

interatomic  behavior  of  a  a  vs  tuts  of  charged  and  neutral  particles  to  the 

macrocopintc  behavior  of  waves  in  interstellar  spaces.  Serious  work  is 

being  done,  based  on  piaern  theory,  to  obtain  a  controlled  and  controllable 

nuclear  fusion  pile.  Plasma  theory  has  influenced  both  semiconductor  and 

gaseous  electronics.  Flaaraa-tnagnef-o-hydrodynamic  based  energy  conversion 

seems  to  ba  within  t  <s  realms  of  large-scale  applicability. 

The  mathematical  formulation  of  plasma  phenomena  in  based  on 

Maxwell’*  laws  of  electromagnetic  field  theory  and  the  laws  of  mechanics. 

Through  a  completely  stochastic  formulation  of  particle  distribution 

functions  in  the  six-dimensional  phase  space  of  mechanics  one  arrives  at 

1  9 

the  Maxwell-Klimontovich  system  ‘  ,  Bv  ensemble-averaging  the  system  and 

assuming  small  fluctuations  one  obtains  thu  non-linear  Haxvell-Vlasov 
1  2 

system  ’  .  Taking  the  moments  of  the  particle  distribution  function  with 
respect  to  velocity  the  fluid  model  is  obtained  which  treats  the  plasma 
like  a  gas  of  charged  particles.  Since  only  a  finite  number  of  moments 
can  be  used  the  fluid  model  cannot  contain  the  same  amount  of  information 
us  the  Uaxwell-Vlasov  or  the  Naxwell-riii.'ontovich  system;  it  restricts 
itself  to  a  few  macroscopic  observables. 


NOT  REPRODUCIBLE 


2 


The  kinetic  model  is  non-linear  as  is  the  fluid  model  in  an 

Eulerian  frame  of  reference.  The  fluid  model  consists  of  a  system  of 

partial. differential. equations  while  the  kinetic  model  is  more  complicated 

and  consists  of  a  system  of  portial-diffcrenciai-intef>ral  equations.  In 

a  Laorangian  reference  frame  the  fluid  model  is  linear;  during  back  crans- 

fo rinat-ion  to  the  Eulerian  frame  however  singular i t ies  prevail  except  under 

3 

the  most  favorable  circumstances  .  One  ia  therefore  obliged  to  look  for 
the  solution  of  ncu-linuar  systems  with  all  the  difficulties  this  presents. 

Since  the  fluid  model  represents  essentially  macroscopic  data,  it 
cannot  explain  phenomena  which  appear  on  a  macroscopic  level  but  can  only 
be  explained  by  the  microscopic  behavior  of  the  system.  Tire  classical  case 
is  the  Landau  damping  which  cannot  be  understood  with  the  fluid  model  alone. 
The  approximations  and  limitations  inherent  to  the  fluid  Model  are  well 

4  5 

known  and  have  bean  well  explained  tv  Vandenplas  and  Gould  and  by  Obertnann  *  . 
Attempt*  to  Improve  the  fluid  model  have  not  been  abandoned;  Feix  proposed 
to  uatt  the  Fourler-Jlarmits  transform  to  obtain  fluid  equations.  A  better 
link  to  kinetic  theory  was  shown,  however  certain  limitations  in  this 
model  fill  exist.  For  example,  such  important  phenomena  as  particle- 
wnve  interaction  cannot  be  studied  with  the  fluid  model. 

The  mathematical  analysis  of  the  kinetic  model  is  far  from  trivial 
oven  when  studied  in  the  linearized  form.  While  for  the  linearized  fluid 
mndui  spectral  analysis  after  Fourier-Laplace  transform  of  Bpace  and  time 
variables  yields  a  finite  number  of  modes  for  a  homogeneous,  isotropic, 
uniform  and  unbounded  plasma,  the  linearized  kinetic  model  reveals  the 


. . "  3 . .  . 

existence  of  a  continuous  spectrum  covering  the  real  frequency  axis 
7 

(Van  Kampen  nodes  )-.  To  e<^lieet-e  matters,  the  real  frequency  axis 
is  the  boundary  of  the  region  of  convergence  of  the  Laplace  transform 
and  Is  therefore  a  branch  of  the  complex  frequency  domain.  Hilbert 
transform-Pleraelj  decomposition  techniques  therefore  are  required  if  one 
wants  to  work  with  real  freouencies.  One  is  confronted  with  the  inherent 
difficulties  of  a  system  which  contains  a  singular  integral  even  after 
the  algebraization  effected  by  the  Fourier-Laplace  transform  of  space 
and  time. 

If  these  are  the  difficulties  which  confront  the  scientist  in 
the  case  of  an  homogeneous,  isotropic,  uniform  and  unbounded  plasma  after 
linearization  of  the  system,  the  difficulties  become  compounded  when  the 
plasma  becomes  bounded  and  is  non-uniform  and  fluctuations  become  large 
enough  to  invalidate  the  linearization.  Various  perturbation  methods 
have  been  used  to  study  the  non-linear  system  when  the  fluctuations  are 

g 

small  ,  Elnaudi  and  Sudan  have  reviewed  available  perturbation  techniques 

Q 

in  their  several  forms  and  pointed  out  their  limitations',  Marcuvitz  has 

utilized  the  two-tine  method  to  studv  plasma  turbulence  phenomena10, 

11  1°  1 3 

Barston  *  “,  Crawford  and  Kino"  have  used  variational  techniques  to 

study  non-uniform  cold  plasma  systems.  Trocheris  has  studied  the  implicit 

conditions  in  the  "ouasi-linoar”  equations  using  both  the  two-time  method 
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and  variational  principles  .  In  general  these  methods  have  only  a  limited 
range  of  validity. 

Most  of  the  work  done  in  the  abov  .  direction  has  been  formulated 
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iii  the  (k,«)  formulae  lor*  (spin-space)  of  the  problem.  The  reasons  for 
this  procedure  are  not  on!''  and  exclusively  the  above  mentioned  mathematical 
fact  that  “transform  "techniques  change  a  differential  system  "into  an  algebraic 
one.  On  the  experimental  side,  due  to  technical  circumstances,  interest 
was  constrained  either  to  oscillatory  phenomena  or  to  aspects  of  long 
duration  with  respect  to  a  plasma  period  T  *»2JT/ai  .  The  results  therefore 

p  p 

hav«  been  mostly  expressed  in  terms  of  spin  space  variables. 

The  investigation  of  strongly  non-linear  phenomena  was  not  only 
hampered  by  the  usual  mathematical  difficulties  posed  by  non-linear  systems. 
These  were  compounded  by  the  fact  that  laboratory  experiments  are  conducted 
in  bounded  and  very  often  non-uniform  plasmas.  The  corresponding  mathematical 
systems  become  extremely  cumbersome  even  in  the  linear  approximation.  Com¬ 
parison  between  experimental  and  theoretical  results  are  often  made  only 
qualitatively.  This  general  difficulty  is  reflected  in. the  literature. 
Comprehensive  treatises  on  plasmas  became  available  since  the  late 
fifties15*16.  Vandonnias  basic  text  on  bounded  plasmas17  lagged  by 
almost  ten  years  and  even  so  restricts  itself  almost  completely  to  the 
fluid  model;  the  kinetic  model  is  used  for  basic  considerations  and  to 
show  that  the  fluid  model  is  a  satisfactory  representation  of  the  system. 

Recently,  ^arcuvitz  was  able  to  make  the  crossover  from  the 
frequency  domain  to  the  tire  domain  from  the  linearized  kinetic  model. 

For  this  purpose,  he  constructed  the  complete  eigenbasis  of  the  Fourier- 
Laplace  transformed  eauationa'10  and  this  in  spite  of  the  fact  that  the 
eigenmodes  are  defined  on  a  branch  of  the  complex  frequency  plane  and  thus 
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need  Plemelj  formulae  for  definition.  On  this  eigenbaee  a  very  convenient 

formulation  for  the  recently  discovered  temporal  plasma  wave  echo  was 

obtained.  The  echo  was  formulated  through  perturbation  expansion  in  such 

a  fashion  that  the  n-th  order  toxin  depends  recursively  or.  the  terms  up  to 

order  fu-1)  but  does  not  contain  the  independent  variable  explicitly  and 

19 

thus  presents  no  secular  terms  .  Unfortunately,  obtaining  even  the 

first-order  term  for  a  Maxwellian  plasma  by  numerical  methods  proved  to 

be  prohibitive  in  terms  of  computer  time  even  on  a  Univac  1103  or  CDC 

6600  ;  the  Fourier  transform  on  a  digital  computer  is  extremely  slow 

and  has  remained  so  in  spite  of  all  efforts  to  the  contrary.  Narcuvitz 

also  attempted  to  obtain  an  elpenbasis  in  the  second  Riemann  sheet  defined 
21 

by  the  problem  ,  The  corresponding,  set  of  eigenvalues  consists  of  enu- 
roerably  infinite  complex  values  with  an  accumulation  point  at  infinity; 
this  set  (whose  least  damped  mode  incidentally  is  the  Landau  mode)  has 
not  yet  been  completely  determined.  Both  eigenbases  were  proven  to  be 
complete. 

The  advent  of  the  high-speed  '-’ipital  computer  changed  the  sit¬ 
uation  to  some  extent.  Not  only  could  one  attempt  the  solution  of  non-linear 
systems,  be  they  algebraic,  differential,  integral  or  mixed  types  by 
numerical  methods;  it  also  became  possible  to  associate  numerical  tech¬ 
niques  with  analytical  results  and  thus  investigate  at  least  numerical 
behavior  in  the  range  of  validity  and  to  point  to  directions  for  further 
theoretical  and  experimental  investigations.  In  the  earlv  sixties,  Dawson 
developed  a  model  in  which  the  motion  of  charged  "sheets"  in  t  slab  bounded 


by  two  parallel  walls  Is  follow  ad;  The  movement  of  the  sheets  is  governed 

2?  £ 
b*/  the  basic  laws  of  physics  .  This  model  was  analyzed  in  depth  by  Feix 

23 

and  has  continued  to  be  used  by  others  ,  Results  proved  to  be  in  sat¬ 
isfactory  agreement  with  experimental  data  in  the  linear  region.  It  is 
presently  being  further  investigated  by  Brooks  and  Cheo  of  the  electro- 
physics  group  at  New  York  University  on  a  nonuniform  slab  and  it  is  hoped 
that  results  can  be  extended  into  the  strongly  non-linear  region  . 

Several  attempts  were  made  to  obtain  solutions  to  plasma  problems 
in  the  time  domain.  Armstrong*^  studied  the  relaxation  problem  of  a  non- 
uniform  electron  plasms  with  a  uniform  ion  background  in  a  plasma  slab 
(see  Equations  1.8  and  1.9,  Appendix  I)  using  Fourier  and  Fourier-Hemite 
series  for  space  and  velocity  respectively  since  these  yielded  the  longest 
time  span  before  numerical  instabilities  began  to  prevail.  As  expected, 
the  smaller  the  deviation  of  the  initial  electron  distribution  the  better 
the  solution  since  the  small  disturbance  parameter  is  needed  to  obtain 
rapid  convergence  of  the  Fourier  series.  A  stable  scheme  for  the  problem 
was  not  obtained  but  results  agreed  with  analytical  data  in  the  linear 
and  quasi-linear  region.  The  solution  for  strongly  non-linear  problems 
was  not  adduced. 

During  1970,  a  very  high  amplitude  and  very  fast  pulse  generator, 

called  bouncing  ball  generator  (SBC) ,  became  available  at  New  York  University 

This  generator  was  constructed  by  P.  Pleshko  in  connection  with  his  Ph.D. 

»  pg 

thesis  work  supervised  ly  Prof.  I.  Palocz  at  New  York  University  »  The 
BBC  is  capable  of  generating  pulses  having,  3db  width  of  approximately 
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1.2  x  10~^  sec.  with  an  amplitude  of  well  over  10^  volts  (appro,  1.5  KV) 
into  a  50  ohm  load  (peak  power  appr.  45  KV7)*  This  capability  of  combined 
high  amplitude  'and  short  duration  and  the  use  of  high  resolution  sampling 
scopes  have  enabled  a  aeries  of  transient  studios  on  laboratory  plasmas 

at  both  linear  and  ncn-llnear  levels.  One  of  the  structures  chosen  for _ 

the  experimental  work  is  the  Tonks-Dattner  structure.  In  this  structure 
a  plasma  column  was  placed  between  a  parallel  plate  system  which  is  con¬ 
nected  to  the  BBG.  The  pulses  from  the  BBG  (width  appr.  0.1  of  a  plasma 
period)  impart  ijfficient  energy  to  the  electrons  and  result  in  stimulated 
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oscillations  of  velocities  reaching  10  m/sec,  well  into  the  non-linear 
region.  Such  a  system  has  great  advantages  over  microwave  C5?  excitations 
since  the  average  power  is  very  small  inspite  of  the  high  amplitude  oscil¬ 
lation.  Thus  the  problem  of  additonal  heating  is  avoided.  Furthermore, 

the  observations  are  made  within  a  short  period  following  each  pulse 
-8 

(appr,  5.0  x  10  sec-.).  This  is  short  comparod  with  average  collision 

times  of  the  system.  Thus  no  ionisation  or  recombination  can  cske  place. 

This  situation  offers  an  attractive  situation  for  theoretical  study;  The 
output  of  the  BBG  can  be  reasonably  wall  simulate  oy  a  Dirac  delta-function; 
colllslonal,  ionization  and  recombination  nodela  are  not  necessary  in  this 
context. 

It  is  the  purpose  of  this  work  to  present  the  results  of  a  theoretical 
effort  which  to  some  degree  parallels  experimental  developments  obtained  and 
being  obtained  at  the  Electrical  Engineering  Department  of  Mew  York  University. 


For  this  stud  •,  the  fallowing  problem  is  considered:  A  non- 
unifortB  kinetically  described  plasma  in  equilibrium  la  bounded  by  two 
parallel  planes  and  contained  by  an  imposed  electric  field.  The  plasma 
is  excited  by  an  impulsive  electric  field  (in  the  fora  of  a  Dirac  delta- 
function)  uniform  in  space  and  perpendicular  to  the  planes.  The  resulting 
phenomena  and  responses  are  obtained  by  solving  numerically  the  system 
formed  by  Maxwell’s  divergence  law  and  Boltzoann-Vlaaov  equations.  The 
necessary  calculations  have  been  performed  on  a  high-speed  digital 
computer  using  a  numerical  scheme  which  has  been  developed  for  this  pur¬ 
pose.  Results,  in  particular  the  electric  field  and  voltage,  are  studied 
in  the  time  domain.  Subsequent  Fourier  transforms  lead  to  frequency 
domain  data  which  are  than  analyzed  and  discussed.  Of  particular  interest 
are  Tonks-Pattner  resonances  in  the  voltage  spectrum  which  are  shown  to 
originate  in  definite  rerions  of  the  slab. 

In  the  next  section,  the  model  is  defined,  formulated,  analyzed 
and  discussed  from  the  mathematical  point  of  view.  Boundary  conditions 
and  initial  conditions  are  established.  Section  III  presents  the  develop¬ 
ment  of  the  numerical  scheme  used  for  the  necessarv  computations  and  dis¬ 
cusses  this  numerical  scheme.  Section  IV  presents  obtained  results  for 
relatively  small  impulses  resulting  in  responses  which  can  be  considered 
linear  to  a  good  degree  of  approximation.  Also  discussed  is  the  effect 
of  variation  of  parameters  necessary  to  define  the  plasma.  Finally,  in 
Section  V,  results  are  adduced  for  large  impulse-excitations  and  some 
observed  non-lir.earities  are  discussed. 
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II.  MATHEMATICAL  ANALYSIS 


A.  General  Coftalderafciong 

From  the  mathematical  point  of  view,  one  is  confronted  in  the 

general  ease  of  a  plasma  con  teining  _n  kinds .  of.  P-8El;  Icles-  with  a-aygtgg.-af - 

S+n  scalar  equations,  four  of  which  are  partial-differentlal-lntegral  of  the 
first  order,  the  others  being  partial  differential  of  first  order.  These 
equations  relate  6+n  dependent  scalar  variables  to  seven  independent  variables 
(See  Appendix  I,),  The  problems  of  existence,  uniqueness  and  well-posedness 
are  tremendous  and  in  order  to  attempt  a  solution  at  all,  even  by  numerical 
methods,  one  has  to  reduce  the  problem  to  simpler  proportions.  One  considers 
therefore  a  system  with  2  kinds  of  particles  only,  electrons  and  ions.  One 
assumes  further  that  the  time  span  considered  in  the  problem  is  relatively 
short  and  therefore  the  ions  can  be  considered  immobile. 

While  it  would  be  desirable  to  obtain  a  solution  in  cylindrical  co¬ 
ordinates  since  most  experimental  work  is  done  in  this  geometry,  it  is  the 
plasma  slab  which  offers  the  minimum  significant  number  of  both  dependent 
and  independent  variables.  The  basic  equations  of  the  kinetic  model  for  the 
pragma  slab  have  been  derived  in  Appendix  I.  It  should  be  noted  that 
Equation  1.9  represents  the  non-linear  Boltzmann -Vlasov  equation  while 
Equation  1.8  represents  Maxwell's  divergence  law.  The  electrostatic  approxi¬ 
mation  is  not  necessary  to  the  system  1.8,9. 

The  first  question  is  what  should  be  the  conditions  at  the  boundary 
of  the  slab?  In  the  experimental  work  the  plasma  is  contained  by  physical 
walls  together  with  the  plasma  sheath;  an  appropriate  simulation  of  this  is 
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given  by  an  externally  imposed  electric  field  (associated  with  an  ion  distri¬ 
bution)  which  rises  rapidly  to  infinity  at  the  walls.  This  high  field  in¬ 
tensity  extends  over  only  a  very  short  distance  (of  the  order  of  less  than 
one  Debye  length)  from  the  wall.  This  field  therefore  simulates  both  the 
impenetrable  wails  and  the  sheath.  (For  ease  of  language,  the  parallel  planes 
which  are  the  boundaries  of  the  plasma  will  ha  called  "walls" although  their 
physical  existence  is  immaterial  to  the  present  model  due  to  the  presence 
of  the  electric  field.) 

The  initial  conditions  are  stipulated  as  follows:  Suppose  that  at 
t»Q~,  the  electron  distribution  function  is  given  by  f(x,v,0-).  At  time 
t^O  one  applies  an  impulse  electric  field  in  the  form  of  a  Dirac  delta- 
function.  Assuming  a  non-relativlstic  formulation,  this  impulse  will  accrue 
the  velocity  of  every  electron  by  an  increment  of  u  in  the  negative  direc¬ 
tion.  In  other  words,  if  an  electron  at  f»0~  had  a  velocity  v,  at  t*(H-, 
this  electron  will  have  a  velocity  v+u.  Therefore  the  initial  condition  of 
f(x,v,t)  at  t“0+  ie  f(x,v+u,Q-).  f(x,v,0~)  is  considered  an  equilibrium 

state  of  the  system. 

The  electric  field  vector  E(x,t)  depends  on  the  charge  distributions 
through  Maxwell's  divergence  equation.  Since  no  particle  has  time  to  move 
from  t*0-  to  t«0+,  the  electric  field  is  continuous  at  t*0. 

To  summarize,  one  considers  the  following  Maxwell-Vlasov  system  as 
a  mixed  initial  value  -  boundary  condition  problem  (formulated  in  MKS 
International  units) 
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defined  on  tbs  domain 
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subject  to  the  initial  conditions 
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and  the  boundary  conditions; 
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The  functions  E*,  f*  and  n*  of  (1)  are  expressed  in  the  MKS 
International  (or  a  similarly  rationalized)  system.  It  is  more  convenient 
to  normalize  by  introducing 
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Kq  is  a  particle  number  density  and  thus  is  a  plasma  angular  velocity. 
(Nq  will  be  chosen  to  be  the  electron  density  at  x*0  in  an  equilibrium 
state  and  u^o  thus  will  be  the  plasma  angular  velocity  at  jp-0,  t«0-  of 
electrons . ) 

!> 

The  problem  is  of  the  mixed  type  in  the  sense  that  both  boundary 
conditions  and  initial  values  are  specified.  They  are  still  expressed 
by  (2)  but  without  the  asterisk.  Physically  expressed,  the  boundary  con- 
ditlons  specify  that  at  no  time  there  are  electrons  on  the  walls  of  the 
slab  nor  do  electrons  with  infinite  velocity  exist. 

It  is  possible  to  gain  some  insight  into  some  properties  of  the 
system  (7,8).  A  short  analysis  has  been  performed  and  it  can  be  shown  that 
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f(XjV,e)  can  be  seen  to  be  hyperbololdal  in  the  sense  that  thare  are  two 
characteristic  manifolds  (in  this  case  two-dimensional  planes)  through  each 
point  of  (x,v,t)-space  but  degenerate  in  the  sense  that  the  quadratic 
form  defined  by  the  respective  partial  differential  equation  degenerates 
into  two  straight  linos  and  does  not  define  a  true  conic.  It  can  also  be 
seen  that  any  place  perpendicular  to  the  x-axis  ie  a  characteristic  plans. 
Thus  any  plane  t*constuit  represents,  a  characteristic  plane.  (It  follows 
that  an  initial  value  problem  defined  at  t »0  with  boundary  conditions  de¬ 
fined  as  above  '‘•institutes  a  variation  of  "Goursat's  problem"  which  solves 
initial  value  problems  defined  on  a  characteristic  at  t«0.)  Therefore  the 
initial  conditions  to  be  stipulated  on  a  characteristic  will  have  to  be 
stipulated  with  great  care  since  they  will  constitute  another  constraint 
on  the  problem.  In  fact,  one  has  to  ask  if  the  initial  values  together 
with  the  boundary  conditions  in  (2)  do  not  overspecify  the  problem  or  whet 
can  be  done  to  avoid  this  situation.  This  again  might  constitute  another 
constraint. 

As  mentioned  before,  the  initial  conditions  at  time  t-0+  result 
from  an  equilibrium  state  at  t-0-.  The  investigation  into  this  problem 
is  not  trivial  slnce  one  has  to  expect  that  the  system  (7,8)  will  impose 
conditions  on  this  initial  state  which  is  defined  on  a  characteristic  of 
the  system. 


To  look  for  a  steady  state  of  the  system  (7,8)  one  sets 
3f(x,v,t)/3t  -  0.  (8)  then  yields 
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Xt  is  shown  (See  Appendix  XX)  that  f(x*v)  i«  separable  under  relatively 
loose  restrictions  on  E(x){  i.e.,  that 

f  (x,v)  -  tp<xW(v)  (12) 


and  it  is  shown  ibidem  that  with  K  and  A  being  constants 
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Thus  the  only  distribution  in  the  equilibrium  state  is  the  Maxwellian 
(13).  For  convenience  one  stipulates 
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and  therefore  (13)  becomes 
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a  being  the  termal  velocity  of  the  electrons.  Appendix  IX  also  proves 
that  <p(x)  >  0  for  all  x  except  on  the  boundary. 

One  further  choses  'RO)  *  1.  Using  (3,12)  the  unnormalised 
equilibrium  electron  distribution  function  is  f*(x,v)  ■  NQV(x)^(v)  and  Nq 
can  be  seen  to  bo  the  number  of  electrons  per  unit  volume  at  x»0  in  the 
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equilibrium  state.  Therefore  defined  by  (6)  has  indeed  become  the 
plasma  angular  velocity  at  x»0  in  the  equilibrium  state. 

Appendix  H  provas^that  fbr  fll^true  one  has,  under  suitable 
mathematical  conditions 
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with  Aqo  the  Debye  length  at  the  origin.  Since  the  electron  density 
varies  with  x,  so  will  the  Debye  length  A^(x).  In  particular,  since 
f(x,v,t)«  0  on  the  walls  of  the  slab,  Ap(x)  becomes  infinity.  This  varia¬ 
tion  of  the  Debye  length  becomes  significant  from  the  physical  point  of 
view.  The  phenomena  described  by  (7,8)  are  macroscopic  in  the  sense  that 
they  refer  to  behavior  outside  of  a  Debye  sphere  centered  on  a  given  particle; 
thus  one  might  neglect  to  a  certain  point  occurrences  within  one  Debye 
length  from  the  boundary.  If  as  can  be  shown,  n(x)  has  a  local  maximum  at 
the  boundary,  then  within  the  slab  there  exists  at  least  one  point  x^o 
where  n(xfflift)  is  a  local  minimum.  If  one  of  these  minima  occur  near  the 
boundary  of  the  slab  in  the  sense  that  the  distance  y  ■  d-*in^n  be  of  the 
order  of  or  less,  it  is  reasonable  to  study  the  phenomena  within  the 
slab  located  between  +xm^n  end  The  region  between  the  minima  and  the 

walls  can  be  assumed  to  represent  the  sheath. 

The  voltage  V(x)  follows  from  (16)  by  integration.  If  one  assumes 
for  the  arbitrary  integration  constant  V(G)  «  0  one  has 
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V(x)  •  -  a 2 An  <p  (18) 

X£  9<x)  is  even  and  moiutonically  depress Jr-g  from  the  center  to  the  walls* 
one  sees  that  the  voltage  is  even  and  monotonically  Increasing  from  the 
center  to  the  boundary  and,  as  follows  from  the  boundary  condition,  becomes 
infinity  at  the  boundary.  The  electric  field  for  its  part  is  odd  sad  fch-ia 
zero  in  the  center  of  the  slab.  The  behavior  of  E(x)  and  V(x)  near  the 
boundary  cannot  be  predicted  without  more  detailed  information  about  «dx> • 

The  question  arises  which  of  the  possible  equilibrium  states  can  be 
initial  states  which  will  follow  the  prescribed  boundary  conditions  after 
the  perturbing  electric  field  has  been  applied.  As  seen  before,  the  initial 
conditions  have  to  be  prescribed  on  a  characteristic  surface  of  the  system 
and  thus  act  as  a  constraint  on  the  problem,  Physically  speaking,  the 
initial  electric  field  (which  is  the  electric  field  of  the  equilibrium 
state)  has  to  be  such  that  even  with  further  developments  it  is  strong 
enough  to  avoid  contact  of  electrons  with  the  wall  surfaces.  It  has  been 
shown  that  the  electric  field  E(x)  of  an  equilibrium  state  which  is  an  ad¬ 
missible  Initial  condition  has  to  behave  near  the  boundary  with  an  order 
l4e  in  1/y  where  y  is  the  distance  fro.”  the  boundary  and  s  is  positive 
finite.  In  fact,  <p(x)  has  to  be  of  a  given  functional  form.  The  following 
outlines  the  procedure  to  obtain  the  obtained  result: 

Assume  now  that  one  has  a  equilibrium  state  background  for  t  <  0 


given  by  (10)  with  (14) 


f (x,v) 


and  obtains 
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3aE 


no 


is 

dx 


(25) 


integration  with  respect  to  x  yielding 

■|§  •  -  «2  u  9  +  constant  (26) 

ot  po 

Examination  of  (7)  row  permits  one  to  assume  that  also  E(x,t) 
should  remain  time-invariant  on  the  walls  and  all  its  partial  time  deriva¬ 
tives  equal  to  sero.  This  is  based  on  the  fact  that  at  no  time  electrons 
of  whatever  speed  can  reach  the  wall.  Thus  since  tP(x)  *  0  on  the  wall, 
the  constant  in  (26)  vanishes  and  (26)  becomes 


at 


Ct)2  u  <p 
po 


(27) 


Similarly,  one  obtains  the  second  partial  time  derivatives  of  f<x,v) 
and  E(x).  The  recursive  formulae  for  the  higher  order  partial  time  deri¬ 
vatives  of  these  two  equilibrium  state  functions  have  been  also  obtained 
and  it  can  be  shown  that  they  are  of  the  form 


3tn+i 


U 


n 

E 

j-0 


Q  .vJ 
aj 


(28) 


where 


flh.  m  _  R 

9tn+1  n 


w  ■  v  +  u 


5»j  ' 


J29) 

(30) 


(31) 
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R  -  R(u,w*  {<?}„,{£}„)  (32) 

n  n  pc  n  n 


a 

t 

(% 

d<P 

dx»  *** 

dn9, 

**  ’  j  *» 

dx 

(33) 

<*>.- 

(E, 

OE 

5t»  •'* 

3nE  . 

(34) 

From  the  form  of  the  terms  Q  ,  and  R  it  has  been  deduced  that  the  functions 

n  j  n 

cp(x)  which  fulfill  the  requirement  that  all  partial  time  derivatives  of 
the  electric  field  and  the  distribution  function  at  time  t»C+  are  zero 


have  to  be  of  the  form 

<P(x)  •  K  exp(-l/h(x))  K  Constant  (35) 

where  h(x)  has  the  properties 

a.  h{x)  is  even  (36a) 

b.  h(x)  is  positive  everywhere  in  the  slab  and  vanishes 

on  the  boundary  and  on  the  boundary  only  (36b) 

c.  h(x)  belongs  to  the  class  of  functions  with  an  infinite 

number  of  derivatives  (36c> 

-3 

d.  h(x)  goes  to  zero  as  y  ,  s  positive  finite  y»d-x  being 

the  distance  to  the  wall.*  (36d) 

From  (16)  it  can  be  seen  that  E(x)  will  go  to  infinity  as  y’<'s+1^. 

K  is  defined  by  9(0)  -  1  as 

K  *  exp(l/h(0) )  (37) 


*Complete  proof  of  the  development  from  Equation  (21)  to  (36)  as  well  as 
explicit  recursion  formulae  (28-34)  can  be  obtained  from  the  author. 
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As  an  of  fur.criv.-.s  rintiafying  the  above  one  has 

cos3  S  <a*  -  xa)S  s  >  0  (38) 

s  being  finite.  These  functions  have  a  aero  of  order  s  on  the  wall  and 
satisfy  all  conditions  (36). 

A  particularly  attractive  h(>t)  out  o£  (33)  is  of  course 
h(x)  *  a4-**  since  all  its  derivatives  of  higher  than  second  order  vanish. 
However*  preliminary  calculations  show  that  due  to  the  form  of  the  terms 
of  the  derivatives  of  the  corresponding  'f'(x)  which  contain  a  factor  of  the 
form 

<?(x)/hn(x)  (39) 

with  n  smaller  or  equal  to  the  order  of  the  derivative  of  9(x),  the  deri¬ 
vatives,  although  zero  on  the  boundary,  have  a  maximum  near  the  boundary 
which  increases  with  increasing  order  and  which  also  gets  nearer  to  the 
boundary  with  Increasing  order.  This  behavior  unfortunately  affects  the 
numerical  solution.  It  proved  therefore  to  be  of  advantage  to  use  an 
equilibrium  state  where  only  the  first  time-derivative  of  f(x,v,t) 
vanishes  at  t*0+.  This  leads  to  loss  of  electrons  as  time  goes  on.  This 
loss,  however,  has  proved  to  be  negligibly  small  compared  with  other  numeri¬ 


cal  inaccuracies. 
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III.  NUMERICAL  SCHEME 


A.  Development  of  the  Numerical  Scheme 

In  order  to  solve  (7,8)  over  n  reasonable  time  compared  with  s 
plasnra  period  T^,  a  numerical  scheme  is  developed.  One  has  the  truncated 
Taylor  series  (note  that  v  not  w  is  used) 

K  aJ  .  , 

£(x,v,t+At)  »  £(x,v,t)  +  2  ■—[£ (x,v,t)](Atr  ~r  (40) 

J*1  'd? 


K  aj 


E(x,t+At)  £  E(x,t)  +  Z  x£E(sc*t> ] (j£st>^  pr 

J-l  3t3  J* 


for  fixed  x  and  v.  One  uses  K«i  for  the  first  approximation. 

Time  discretization  is  based  on  the  constant  Interval  h»At.  Then 
for  a  non-negative  integer  j 


f(x,v,(j+l]h)  “  f (x*v, jh)  +  h  (x.v.jh)] 


E(x,[j+i]h)  5  E(x,jh)  +  h  fr[E(x,jh)) 

ot 


the  reminders  being  0(h2).  Similarly  x  and  v  are  discretized.  One  has 


therefore 


t  *  Jh 

j  -  0,1,2 . 

...  *«J 

(44a) 

x  ■  mAx 

m  *  0,+l,-t2. . . . 

. .  »+M 

(44b) 

v  -  nAv 

n  ■  0,+l,+2. . . . 

. .  .+ri 

(44c) 

1 — rxm-HsMsvstasHi 


\ 


where  Jh  ia  the  total  time  to  be  conside  ed  and  HAx  *  4,  4  being  the  half¬ 
width  of  the  slab.  H  ia  chosen  to  cover  a  deviation  of 

v  »  4»'*?a  (45) 

n 

from  the  average  velocity.  For  a  GauBsian  this  truncation  introduces 
-  -26. 

an  error  or  0(e  ). 

Goe  now  nets 


f (raix.nAv, jh)  £  tJ 
s  J  m,n 


jh)  5 


and  (40-41)  become 


J+i  „  fi  +  h  f J 

m,n  in ,  n.  bt  m,n 


■'41  i  3  i 

EJ  *■  nJ  +•  h  - —  EJ 

m  m  3t  m 


One  now  uses  (8)  and  sets  using  central  differences 


i_  fJ 

3t  m 


f  J  a  -  — i_  v  ff3  J  )  4.  .  rl  /f.1 

m,n  2ux  n  w-l,n  2Av  m. 


fJ  \ 


n+1  "m.n-l 


One  sets  further 


and  (49)  becomes 


f3  =  D*f^ 
ot  m,n  in,n 


f3+1  -  f3  4-  hD'fj 
m,n  tn,n  m,n 


In  order  to  use  (48),  one  sets  (31)  in  a  discrete  form  with  respect  to 
v  (which  replaces  w) 


) 


V.-  - 


t 


% 


i- 

I 


i 


I 


9ZE 
3x3 1 


4  +N  1 
~  w2  Av  En  D-fJ 


po 


-N 


m,n 


(52) 


and  then  replaces  the  Integral  over  x  by  the  sum 


Vp  j  *  -  w2  AxAv 


t!  +N 


..  J 

v  -  U 


A 


po 


■  (M-l) 


Zn  U-f; 

-/  k‘n 


(53) 


or  alter natively  by  the  summation 


a  .  m  +N  . 

vrJE3!  -  +  w*  AxAv  E.  £„  D-ff 

at  tu  po  k  n  k,n 

M-l  -N 


(54) 


where  the  constant  of  integration  was  set  to  zero  due  to  the  argument 
previously  developed  (See  after  Equation  26).  (53)  and  (54)  show  quite 

clearly  that  for  an  arbitrary  initial  condition  the  imposed  boundary 
conditions  can  overspecify;  (53,54)  might  easily  lead  to  different  re¬ 
sults  for  the  same  m.  If  however  a  proper  initial  condition  is  used,  the 
two  equations  can  be  used  to  reduce  numerical  error. 

If  one  sets  similarly  to  (50) 


(48)  can  be  written 


3_ 

3t 


E 


i 


m 


D-E^ 

m 


(55) 


E^+i  »  eU  h  D-E^  (56) 

tn  m  m 

(51,56)  constitute  a  'Vmeryless"  advancing  in  time  from  jh  to  (j+l)h 
using  exclusively  the  information  contained  in  the  system  at  time  jh. 
From  the  numerical  point  of  view,  however,  "leapfrogging"  seems  to  yield 


24 


better  results  und  better  stability  in  the  sense  to  be  explained  in  the 
next  section.  One  replaces  < 5 1  * 5 ^ >  by 


-■l-l 

1*  -* 

”m,n 


+  21'.  0 


*3 

"n,n 


(5?) 


Kyyl  -  2'j~l  +  2h  D-Ej 

si  m.  m 


(58) 


except  for  j^O  -.Men  (51,56)  continue  to  be  uced. 


E.  Cons! derations  on  N’.sr/»r J  cal  Stability 

There  is  ao  Vocwn  technique  applicable  to  tfyj.  present  problem 
for  a  study  of  stability  criteria  In  the  -snse  of  Vcn  Neumann.  It  is 
possible,  by  msr.iptili.ii ag  tha  system  (7,8),  fcc  arrive  at  e.  quasi- 
linear  p?.rt.iril  differential  equation  of  third  order  for  tha  function 
£ (x , v , l )  vhich  can  be  analyzed  for  Von  Neumann  stability  in  a  heuristic 
sense.  However,  besides  losing  two  levels  of  Information,  the  stability 
analysis  will  only  yield  results  in  a  "local"  sense  offering  no  assurance 
of  global  stability. 

A  large  number  however  of  computer  runs  has  shown  that  a  macro¬ 
scopic  parameter  exists  which  will  define  the  area  in  which  obtained 
results  are  reliable.  This  parameter  is  the  total  number  of  particles 
per  unit  area  of  slab.  This  number  is 


not  reproducible 
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■KJ  +00 

K^Ct)  ■  t?c  /  iii  /  dv  f (x,v,i;) 
-d  -oo 


M-l 

*  No4xAv  lm 
-M+l 


f 


J 

m,n 


“  NT(jh)  (59) 

One  defines  the  "normalised"  electron  quantity 

•  NT(jh)/NT(0)  (60) 

which  ideally  should  be  constant  and  equal  to  one.  However,  besides 
numerical  instability,  there  Is  a  factor  to  be  considered  which  will 
allow  to  vary  at  least  to  some  degree.  This  is  the  replacement  of 
a  continuous  electric  field  £(x,c)  by  the  discretized  field  E(m,t). 

E(x, t)  should  go  to  infinity  with  an  order  (1+s)  as  x  goes  to  the  wall. 
(s»l  was  considered.)  This  Is  impossible  to  simulate  on  any  computer. 

Even  working  with  the  finest  possible  mesh  in  x-space,  the  electric  field 
near  the  boundary  cannot  be  made  infinite.  Therefore  some  electrons  are 
lost  to  the  boundary.  This  is,  in  fact,  observed  during  a  very  small  time 
compared  to  Tp  but  then  begins  to  Increase  (and  sometimes  to  oscillate 
very  slightly).  Computer  runs  were  made  with  several  values  of  M  and  N 
using  either  (51,56)  or  (57,58)  or  any  combination  thereof.  The  physical 
parameters  of  the  plasma  were 
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Nq  »  10“'  electrons  per  cu.  meter  (61a) 

a  •  10^  m/sec  (61b) 

d  -  5.10"3m  (61c) 

<p(x)  *  cos2  (61d) 


and  for  several  values  of  h  ranging  from  T  /20  to  T  /240  and  (where 

P  P 

T  is  the  plasma  period  of  the  center  of  the  slab  in  equilibrium)  for 
F 

several  values  of  u  ranging  from  ,1a  to  a.  As  long  as  was  leas  than 
1.004,  results  for  the  same  parameters  were  practically  the  same  inde¬ 
pendent  of  the  used  scheme.  However  scheme  (57,58)  (which  in  some  way 
is  similar  to  the  Duf ort-Frankel  scheme  of  the  diffusion  equation)  always 
reached  1.004  after  much  more  time  Chan  any  other  scheme,  and  thus  can  be 
taken  as  "more  stable".  After  several  dozen  of  comparative  runs,  all 
schemes  were  discarded  in  favor  of  (57,58).  The  macroscopic  parameter 
evaluated  in  this  stage,  except  N*^,  was  the  voltage  both  wall-to-wall  and 

from  -x  .  to  x  .  . 
min  min 

The  problem  thua  became  to  extend  the  time  (in  terms  of  T  )  at 

P 

which  reached  the  value  1.004.  Call  this  time  N^c.  This  time  is  in 

strong  dependence  of  the  strength  of  the  impulse  given  by  u/a  as  can  be 

expected  since  this  value  is  a  measure  of  the  non-linearity  Introduced 

into  the  system.  For  decreasing  u/a,  NJ  increases  with  constant  h,M,N. 

UJ  also  increases  with  increasing  M  and  U  for  constant  h  provided  h  be 

i  ic 

sufficiently  small.  The  total  amount  of  time  at  which  NJ  ■  also 
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lncresses  for  decreasing  h  with  constant  M  and  N  for  some  time  after 
vhich  a  decrease  of  h  yields  no  further  improvement.  Thus  the  technique 

A  £ 

to  extend  the  time  at  vhich  NJ  is  reached  seems  to  consist  in  increasing 
M  end  N  (interestingly  best  results  correspond  to  M*N)  end  then  decreasing 
h  until  no  further  improvement  is  obtained.  Unhappily,  this  also  requires 
an  increase  in  computer  core  space  and  processor  time.  If  one  doubles  M 
ana  H  and  halves  h,  the  computer  core  requirement  increases  by  a  factor 
of  4  and  CPU  time  by  8.  Although  New  York  University's  Heights  Univac 
1108  is  a  very  fast  and  large  machine,  the  program  is  operating  near  the 
core  capacity  of  the  1108  and  uses  a  considerable  amount  of  CFU  time 
per  run. 

The  increase  of  the  number  of  electrons  was  investigated.  Essen¬ 
tially,  it  consists  in  the  fact  that  in  the  "corners"  of  the  (x.v)-space 
the  numerical  values  of  the  central  differences  in  (49)  show  a  relatively 
large  error  compared  to  true  analytical  values.  The  two  errors  can  rein¬ 
force  each  other,  deform  D‘f^  and  thus  distort  and  then  E-^.  With 

m,n  m,n  m 

time,  these  errors  propagate  through  the  whole  slab  and  lead  to  a  catas¬ 
trophic  explosion  of  .  However,  as  long  as  does  not  pass  1.004  by  a 
considerable  amount,  the  resulting  distortion  seems  to  be  negligible;  it 
certainly  does  not  affect  voltages,  electron  density  and  electric  field 
and  discrepancies  in  velocity  space  are  practically  restricted  to  the 
corners  of  (x,v)-space. 

In  order  to  get  an  idea  about  the  variation  of  and  the  time  at 
lc 

which  NJ  is  reached,  Table  I  shows  a  few  results  obtained  for  the  plasma 
(61).  This  plasma  In  equilibrium  has  at  the  center  of  the  blab  a  frequency 
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TABLE  I 


Values  of  at  times  t/T  for  various  sets 

P 

of  values  T  /At,  d/Ax,  4av^/Av  and  u/a 


10 


0.005 


1 

1 


3 


_d_ 

wT 

Ie 

t 

a 

Ax 

Av 

At 

To 

0.1 

50 

50 

320 

1.004 

8.050 

0.2 

60 

60 

400 

1.004 

5.657 

0.5 

60 

60 

400 

1.004 

2,867 

1.0 

60 

60 

400 

1.004 

2,055 

0.1 

so 

50 

320 

1.0073 

10,0 

0.05 

50 

50 

320 

1.C034 

10,0 

0.01 

50 

50 

320 

1.0005 

10.0 

0.01 

50 

50 

320 

1.0009 

15.0 

0.01 

50 

50 

320 

1.001 

20.0 

0.1 

60 

60 

100 

1.0095 

0.180 

s 

£ 
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of  approximate  lGRz  and  a  temperature  of  approximately  30000°K.  The  first 

4  f 

group  of  results  in  Table  I  shows  values  of  NJ  for  various  values  of  the 

involved  parameters,  in  particular  for  values  of  u/a  from  0.1  to  1.0.  The 

second  group  shows  values  of  for  decreasing  values  of  u/a.  The  last 

set  points  out  that  a  very  strong  sensitivity  seems  to  exist  in  terms  of 

the  relationship  between  H,N  and  h.  This  sensitivity  is  in  accordance 

with  similar  phenomena  observed  for  other  systems  of  partial  differential 

equations  solved  by  numerical  difference  methods. 

The  values  shown  in  Table  I  are  pretty  characteristic.  They  vary 

with  the  half-width  d  although  not  in  straight  proportion.  tr  increases 

with  d  while  decreases  with  constant  t/T  or  t/T  increases  with  ccn- 

P  P 

stant  N* ,  all  other  parameters  being  constant.  It  follows  from  Table  I 
that  unhappily  at  the  present  time  it  is  not  possible  to  Investigate  the 
effect  of  very  strong  impulses  resulting  in,  say,  u/a  ■  1  for  relatively 
long  periods  of  time.  Even  so,  some  non-linearities  resulting  from  the 
increase  of  u/a  up  to  a  value  equal  to  unity  have  been  observed  and  will 
be  reported  in  the  respective  section. 

4  £ 

The  quest  for  extension  of  Ir  also  explains  the  'use  of  (61d). 
This  function  proved  to  yield  larger  R^c  than  (35)  with  h(x)  *  a2  -  x2 
although  (61d)  does  not  satisfy  conditions  (35,36).  As  mentioned  before, 
this  should  result  in  some  (additional)  loss  of  electrons  to  the  boundary 
but  the  observed  numerical  results  in  this  direction  have  been  negligible. 
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IV.  RESULTS  IK  TKE  LINEAR  REGION  DUE  TO 
S’. 'ALL  AMPLITUDE  EXCITATIONS 

A.  Basic  Features  of  Voltage  Behavior  and  the  Fourier  Transform 

In  accordance  with  the  arguments  presented  In  section  II. b 
(following  (17)),  the  voltage  henceforth  will  be  understood  to  be  the 
voltage  between  Che  two  minimum  density  points  located  at  •  Ibis 

avoids  the  numerical  difficulties  caused  by  the  singular  field  at  the 
walls.  It  gives  an  accurate  description  of  the  voltage  between  the 
walls  since  the  singular  electric  fields  cancel  out.  In  the  equilibrium 
state,  the  voltage  between  any  two  points  symmetric  with  respect  to  the 
center  of  the  slab  is  always  zero  and  will  remain  so  at  t*0+  due  to  the 
continuity  of  the  electric  field.  Thus  the  voltage  at  t«G  will  always 
be  zero. 

It  is  convenient  to  introduce  a  "standard  experiment".  The  equi¬ 
librium  plasma  for  this  experiment  will  be  the  plasma  slab  defined  by 
(61).  The  excitation  for  the  "standard  experiment"  will  be  supplied  bv 

an  impulse  electric  field  yieldinp  an  electron  velocity  increment  of 

5  5 

u*10  towards  the  wall  at  x=-d.  The  stimulus  u»lG  has  been  chosen  for 

the  "standard  experiment"  since  this  value  proves  to  be  approximately  the 

upper  limit  of  the  linear  region. 

Some  data  for  the  "standard"  plasma  slab  in  equilibrium  are 

listed  in  Table  2. 

All  numerical  experiments  henceforth  will  be  Performed  on  plasmas 
having  No  and cp (x)  given  by  (61a,d).  The  equilibrium  thermal  velocity  a. 
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TABLE  2 


Equilibrium  Plasma  Data 


Standard  Experiment 


Electron  density  N0  at  center  of  slab: 
Thermal  velocity  a  : 
Electron  density  law  : 
Velocity  distribution  : 


1016  el/m3 


10  m/sec 
N 


„  2  It  x 
cos  7  3 


Csussian 


At  center  of  slab: 

Plasma  angular  velocity  w  : 

Plasma  frequency  f  : 

Plasma  period  : 

Debye  length  : 


5.6452.109  rad/sec 
0. 89846. 109  KHz 
1.1130  n  sec 
0. 177142. 10'3  m 


For  a  slab  with  halfwidth  of  0.005m,  at  +xm<n  (minimum  ion  density  points); 


Electron  density,  relative  to  same  at  center: 

Ion  density,  ditto  : 

Ion/electron  ratio  : 

Voltage  : 

-3 

Electron  Debye  length  *  0.6314.10  » 

x  .  *  0.9048.103  ;  0.18d  sl  1.433  local  Debye 


0.07870 

0.1574 

2.000 

14.45  volts 


lengths 
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the  slab  half-width  d  and  the  stimulus  u  fill  !e  varied  according  to  need. 
For  brevity  everv  numerical  experiment  "ill  be  described  b«  tf’ese  three 
parameters. 

For  easy  comparison  of  results  arid  convenience,  one  defines  the 
normalized  voltage  and  time 


V  .  Ji-. 
n  IQ5 


0.005 

d 


V 


(62) 


tn  -  t/Tp  (63) 

henceforth,  the  terms  "voltage"  and  "time"  will  always  mean  the  normalized 
voltage  and  normalized  time  as  defined.  The  sul  cript  "n"  will  be  dropped. 

Fig.  1  shows  the  behavior  of  the  voltage  over  ten  plasma  periods 
for  the  "standard  experiment".  Table  3  shows  the  values  of  the  voltage 
peaks  and  their  time  of  occurrence.  Computer  calculations  were  performed 
using  100  intervals  in  both  x-space  and  v-space  and  a  tine  interval  of 
1/320  plasma  periods  h«l/320,  see  (44)  et  sen.).  Finer  inter¬ 

vals  were  tried  but  did  not  yield  any  improvement  in  the  linear  region. 

Peak  times  shown  in  Table  3  are  therefore  correct  to  0.003. 

Examination  of  Fig,  1  and  Table  3  shove  a  rather  puzzling  behavior. 
One  would  expect  a  reasonable  periodic  oscillation  damped  according  to  some 
reasonable  damping  lav;.  Examination  of  Fig.  1  reveals,  to  mention  a  few 
aspects,  that  the  third  peaks,  both  positive  and  negative,  are  larger  than 
the  respective  second  peaks.  The  fifth  positive  peak  is  larger  than  the  2nd, 
4th  and  5th  and  almost  as  large  as  the  3rd.  Time  spans  between  adjacent 
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peaks,  hencefort1'.  caller  "neak-to-paak  times",  vary  from  0,703  around 
t*«4.o  to  1,869  f  tpnroximateiy  t*9.  One  observes  further  that  between 
the  fourth  and  fifth  negative  peaks  the  voltage  alwavg  re-'alns  negative; 
this  remains  so  even  after  the  dc  bias  in  the  dnta  is  removed.  Examination 
of  the  voltage  between  the  last  negative  and  posit!  a  peaks  of  Pip.  1 
reveals  that  the  first  difference  passes  through  a  minimum  at  approximately 
t=9  while  remaining  always  positive. 

Tills  apparently  erratic  behavior  could  he  suspected  to  be  due  to 
the  numerical  calculations  aud/or  instabilities.  A  numerical  experiment 
was  performed  with  all  parameters  enual  to  those  of  the  standard  experiment 

4 

except  the  stimulus  reduced  to  u«10  ,  1/10  of  the  standard  experiment. 

The  voltage  (after  normalization)  is  identical  to  the  one  observed  in  the 
prior  case  except  for  a  very  small  difference  around  t-10  where  nunericsl 
instabilities  begin  in  the  case  of  the  standard  experiment  (see  Table  1). 

The  voltage  'or  this  experiment  with  reduced  excitation  is  plotted  in  Fly,  2. 
Extended  observation  is  now  possible.  The  same  "erratic"  behavior  contin¬ 
ues, 

A  numerical  experiment  was  also  performed  for  a  plasma  with  a  and 

4 

d  eoual  to  the  standard  experiment  and  u*=5,lfl  .  results  matched  those  of 
Table  3  and  Figs.  1  and  2.  Results  are  not  shown  since  they  contain  no 
new  information. 

Fourier  transformation  of  the  data  displayed  in  Figs,  1  and  2  pre¬ 
sented  significant  results.  It  is  however  necessary  to  look  at  the  trans¬ 
formation  before  attempting  an  interpretation  of  the  results.  As  far  as 


3B 

the  necessary  numerical  integration  itself  is  concerned,  trial  computations 

showed  that  the  use  of  20  points  per  time  unit  (i.e.,  per  plasma  period) 

yielded  a  sufficiently  good  anproximatler.  (320  points  per  plasma  period 

are  available).  A  seriou3  problem,  however,  is  posed  by  the  fact  that 

the  voltage  is  cnly  brown  for  a  relatively  short  time  period  and  not  to 

infinity.  Therefore  the  integration  will  not  yield  the  Fourier  Transform 

of  the  voltage  but  the  Fourier  Transform  of  the  voltage  multiplied  bv  a 

truncation  function  l!(t)  -  L’(t-t  ),  where  t>(t)  is  Ksavyside’s  step  function 

tn 

and  t  is  the  time  limit  up  to  which  the  function  is  known .  The  result  of 
the  transform  thus  will  be  the  convolution  of  the  voltage  (known  to  in¬ 
finity)  with  the  capital  function.  In  practice,  one  has  to  construct  a 
curve  through  the  local  manir.a  of  the  obtained  transform  of  the  truncated 
voltage.  It  is  also  known  that  the  number  of  local  maxima  per  frequency 
unit  is  approximately  the  number  of  tine  units  during  which  the  voltage  is 
known.  In  other  words,  if  the  voltage  is  known  for  10  time  units,  one  can 
expect  10  local  maxima  per  frequency  unit. 

Since  Fourier  transform  of  voltages  known  over  different  time 
spans  will  have  to  be  performed,  it  becomes  necessary  to  renormalize  the 
time  scale.  The  normalisation  adopted  herein  practically  results  in  the 
fact  that  the  DC  value  of  the  Fourier  transform  eouals  the  arithmetic  mean 
of  the  voltage  in  the  time  domain. 

Frequency  normalization  is  imposed  by  (63).  If  fn0*  1/T^  then 

f  =  f/f  =  f-T  (64) 

n  po  a 

not  kepR°°',c,bU 
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is  the  normalized  frequency.  Again  the  term  "frequency"  henceforth 
will  imply  normalized  frequency  and  the  subscript  ' V  again  will  he 
dropped.  {The  angular  velocity  t:  however  will  not  he  normalized).  All 
tables  and  graphs  of  Fourier  transforms  display  the  amplir  " .  vcrr.us  nor¬ 
malized  frequency. 

In  order  to  study  the  effect  of  truncation,  the  Fourier  trans¬ 
forms  using  the  uunet i.cai  equivalent  to  Fig.  2  v’cre  norforncd  based  on 
the  first  10,  the  first  15  nrui  the  first  20  time  units,  results  are 
listed  in  Tables  4  ar.d  shorn  in  Fir.  3.  First  one  notes,  before  attempting 
to  interpret  the  results,  that  the  agreement  between  the  three  curves  in 
excellent;  the  difference  in  magnitudes  results  from  the  normalization 
process  together  with  the  attenuation  to  be  expected  with  time.  The  main 
difference  between  the  transform  based  on  10  tine  units  as  compared  to 
the  other  two  is  that  the  valley  around  f»0.9  does  not  show  up.  This  ic 
due  to  the  scarcitv  of  local  maximum  points.  Examination  cf  the  transform 
in  this  region  shows  the  existence  of  a  minimum  value  of  the  phase  at 
which  the  absolute  value  has  a  magnitude  very  near  to  the  one  show^n  by 
the  other  transforms.  rnenever  the  scarcity  of  local  maxima  will  force 
the  use  of  this  procedure  in  order  to  obtain  points  of  the  Fourier  trans¬ 
form,  this  will  be  noted  in  the  tables  as  "phase  reversal  points"  and 
this  term  also  shall  be  used  in  the  text  whenever  necessary. 

In  Table  5  the  Fourier  transform  envelope  points  are  listed  for 
the  standard  experiment.  Graphical  representation  is  given  in  Fig.  A. 
Agreement  with  data  of  Tabic  4  and  ^ig.  3  is  excellent.  T’ote  in  particu¬ 
lar  the  phase  reversal  point  at  f =0.855, 

NO1?  RtP 


H 


rQDUC\BU 


1*0 


TABLE  l* 


Fourier  Transform  of  V(t) 

,  over  10, 

15  and  20  * 

r 

P 

a  •  10^ 

u  "  10^ 

d  -  .005 

20 

T 

_ E 

.000 

.1797 

.030 

.2014 

.075 

.1997 

.125 

.2174 

.180 

.2157 

10  T 

15 

T 

.235 

.2370 

_ £ 

_E 

.295 

.2362 

.000 

.2855 

.000 

.2006 

.330 

.2536 

.050 

.4872 

.035 

.3103 

.385 

.3007 

.150 

.4918 

.100 

.3222 

.440 

.3742 

.250 

.5450 

.165 

.3031 

.500 

.5951 

.345 

.6351 

.230 

.3352 

.580 

1.4761 

.450 

.8540 

.300 

.3632 

.655 

.1352 

.570 

2.0878 

.360 

.3900 

.735 

.5357 

.695 

.6610 

.425 

.4330 

.800 

1.3965 

.800 

1.6282 

.495 

.7435 

.865 

.7742 

.985 

2.2913 

.580 

1.8835 

.925 

.9057 

1.115 

.8188 

.655 

.3413 

.590 

1.4140 

1.225 

.4801 

.720 

.7155 

1.115 

.5163 

1.330 

.3420 

.800 

1.4657 

1.175 

.3673 

1.430 

.2663 

.875 

.6747 

1.235 

.2686 

1.530 

.2100 

1.000 

1.6245 

1.320 

.1857 

1.630 

.1735 

1.095 

.7668 

1.380 

.1852 

1.730 

.1496 

1.165 

.4090 

1.435 

.1587 

1.330 

.1281 

1.23C 

.3250 

1.520 

.1024 

1  930 

.1147 

1.300 

.2517 

1.580 

.1257 

2.C30 

.1033 

1.365 

.9889 

1.535 

.1091 

2.130 

.0913 

1.435 

.1733 

1.590 

.0652 

2.225 

.0816 

1.500 

.1405 

1.725 

.0748 

2-325 

.0781 

1.565 

.1333 

1.780 

.0913 

2.425 

.0683 

1.635 

.1141 

1.335 

.0826 

1.700 

.0992 

1.890 

.0550 

1.765 

.0891 

1.920 

.0548 

1.830 

.CS11 

1-5  0 

.0733 

1.900 

.0782 

2.055 

0692 

1.965 

.0680 

2.090 

.0456 

2.035 

.0631 

2.125 

.0372 

:.:oo 

.0575 

2.180 

.0577 

2.165 

.0580 

2.235 

.0559 

2.235 

.0514 

2.290 

.0358 

.0472 

2.  '35 

.0351 

2.365 

.0435 

2.700 

.0458 

2.  i  30 

.0443 

2.435 

.0492 

2.500 

.0429 

2.495 

.0331 
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TABIE  5 


FIG.  4  FOURIER  TRANSFORM  OF  V(t) 


uu 


Results  shown  graphically  in  Figs,  1-4  together  with  results 
for  ft  plasma  ala'  with  identical  a  and  d  and  u  indicate  that 

vithln  this  region  the  system  operates  linear lv,  Variation  of  the  re¬ 
sponse  as  a  function  of  a  and  d  fill  be  discussed  in  another  section. 

A  tvplcal  behavior  of  Tonks-Dattner  resonances  is  thus  observed. 
Specifically,  three  high-amplitude  peaks  can  be  seen  at 

1  -  f  s  1.0C0 

2  -  f  -  c.eoo 

3  -  f  S  0.575 

In  addition,  another  set  of  resonances  of  smaller  amplitude  can 
be  seen  to  exist  at  higher  frequencies.  The  voltage  for  this  range  of 
frequencies  Is  displayed  in  Fir.  5.  Those  resonances  clearly  are  thermal 
resonances.  They  do  not  shov  up  in  other  voltages  herein  since  the  number 
of  points  per  unit  of  frequency  is  too  small  to  provide  the  necessary 
detail,  'hlle  these  resonances  are  easily  explained  and  the  resonance 
at  approximate!'--  unit  frequency  is  the  tain  resonance,  the  resonances 
appearing  at  frequencies  of  approximately  C.8  and  C.575  necessitate  fur¬ 
ther  investigation.  In  order  to  present  the  necessarv  clarification,  i  ■ 
is  logical  to  studv  the  behavior  of  the  electric  field  vith in  the  plasma 
slab.  This  is  carried  out  in  the  next  section. 

Some  thermal  resonances  at  frequencies  above  unity  can  be  observed 
in  Fourier  transforms  of  the  electric  field  or  its  deviation  adduced  in 
the  next  .:ection.  Since,  hoi -ever,  they  are  not  the  o'-ject  of  this  study, 
no  further  attention  need  be  given  to  the  sare. 


NOT  REPRODUCIBLE 
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B.  The  Electric  Field  and  its  Fourier  Transform 

In  this  section,  the  behavior  of  the  electric  field  «»ill  be 
studied.  Since  the  equilibrium  electric  field  resultc  in  zero  voltage, 
the  investigation  rill  be  centered  on  the  deviation  AE(x/d,t).  Results 
of  computer  calculations  have  sho*m  that  AE  is  approximately  even.  Thus 
studying  this  function  in  one  half  of  the  slab  rill  be  sufficient.  Data 
obtained  are  shown  and  plotted  oniv  for  -d<x/d<0.0  unless  noted  (the 
si^n,  however,  will  he  omitted  whenever  convenient). 

The  deviation  of  the  electron  density  An  is  also  studied.  This 
deviation  is  proportional  to  -A(3E/3x)  ns  follows  from  (7)  and  yields 
therefore  information  on  the  behavior  of  AE.  Further,  the  electron  migra¬ 
tion  behavior  is  also  of  interest.  The  function  An  is  odd  in  x/d  to  a 
verv  pood  approximation  and  thus  practically  zero  at  the  center  of  the 
slab  in  the  linear  revion.  (In  fact,  one  can  almost  define  the  linear 
repion  by  this  propertv  of  the  electron  density  deviation  at  the  center 
of  the  slab.) 

Fig.  (>  shows  AE  vs.  time  for  x/d«-Q.O,  0.5,  O.S  for  the  standard 
experiment.  Its  three  graphs  bear  no  resemblance  to  each  other.  At  the 
center  of  the  slab,  x/d«0.0,  the  el  ctric  field  is  an  attenuated  sinusoid 
at  a  Ereiucncy  near  unJtv,  i.e,,  the  local  plasma  freruenev.  Further 
towards  the  wall,  at  x/d»0,5,  inltiallv  ponU-to-pcr.h  times  arcs  slightly 
above  unltv,  then  become  about  1.5,  go  back  to  unity,  increase  again. 

Nunr  the  wall,  at  x/ri«0.8,  peak-to-penk  times  are  about  1,5  except  for 
the  ’Vlpple"  between  t»4,0  and  t»5,Q, 
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To  look  at  the  phenomena  from  a  different  viewpoint.  Fig.  7  plots 
the  quantized  values  of  An  vs.  x/d  and  time,  showing  the  electron  nigra— 
tion  pattern.  Twelve  discrete  levels  of  quantization  are  chosen  desig¬ 
nated  by  ME",  "D",  "C",  "B",  "A",  "1",  " 2 ",  "3",  "4" ,  "5"  in 

decreasing  order  with  the  zero  located  between  'V  and  Occasionally 

the  lecter  "X"  is  used  to  designate  a  level  sufficiently  near  zero. 

As  mentioned  before,  An^O  at  x/d»0.0.  One  observes  that  the 
first  and  second  "positive"  and  the  first  "negative"  peaks  travel  froo 
the  center  of  the  slab  towards  the  sheath.  The  second  "negative"  peak, 
however,  docs  not  reach  x/d“0.8  anymore.  It  seems  to  dissipate  Itself 
at  x/d»0.6,  0.7  around  t»2.3.  Some  later  peaks  arc  seen  to  travel  from 
the  plasma  sheath  towards  the  center.  One  notes  on  almost  continuous 
zero  between  x/d-0.5  and  x/d®0,6  from  t*3. 5  to  t»5.5.  Cf  interest  is  also 
the  "bunching"  of  particles  at  x/d«0.5  where  An  remains  positive  from 
t«5.2  to  t*6.4, 

Fip,  «  displays  AE  in  a  similar  fashion.  Similar  to  Ftp,.  1, 
one  can  see  the  first  three  electric  field  "peaks"  (here  the  first  one 
is  "negative")  travel  towards  the  wall.  The  fourth  peak  (second  "posi¬ 
tive")  which  3torts  at  the  center  around  t»l. 75  does  not  reach  the  \*all. 

At  x/d "0.6  there  exists  for  about  0.5  time  units  n  barrier  which  separates 
positive  and  negative  values  of  AE  at  about  t«2.5  ant’  on  both  sides  elec¬ 
tric  field  wnvea  seem  to  develop  independently  of  each  other.  The  posi¬ 
tive  pend;  near  the  wall  then  begins  to  travel  to'ards  tho  center  and  met' 
pas  with  another  potiit  ivo  peek  coming  from  the  center  of  tho  slab.  Another 
"harrier"  can  be  observed  at  x/d»ft,3  from  t»4, 3  to  t«?,0.  The  oscillation 
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FIG.  Bb:  QUANTIZED  &E(t) 

t/T  =0.0(0. 05)7-5,  x/*HLG{0.i50.8 
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of  AE  at  this  point  of  the  slab  is  small,  while  both  towards  the  center 
and  towards  the  wall  more  intense  chances  can  be  observed. 

Fig.  9  presents  in  a  different  way  the  information  contained 
in  Figs,  1  and  8,  AE  and  An  are  plotted  va,  x/d  at  12  different  points 
per  time  unit  for  0.0<u<3.0  and  6.1<t<7.5.  Fig.  9  permits  one  to  study 
the  periodicity  of  AE  and  An  throughout  the  slab  and  it  becomes  easy  to 
sea  the  differences  of  behavior  in  the  two  shown  intervals.  For  Instance, 
one  note.s  that  the  oncillationo  of  An  in  the  region  from  x/d»G.C  to  x/d»0.3 
have  become  small  with  respect  to  the  ones  in  the  rest  of  the  slab  during 
6.!<t<7.5,  Fig.  9  also  permits  one  to  follow  the  travel  of  the  peaks  in 
a  clearer  fashion  than  Figs.  7  and  8  especially  the  quantitative  cspccts 
involved  in  this  "travel". 

Analysis  o’  Figs.  6-f  points  to  the  existence  of  three  character** 
is  tic  regions  within  the  olab  (or  five  separata  sectors,  t»?o  pairs  of 
which  are  symmetric  with  respect  to  the  center  of  the  slab,  each  pair 
forming  a  region). 

This  is  most  dramatically  demonstrated  by  the  Fourier  arectra 
of  the  electric  field  (dynamic  part)  shown  in  Tnl  lo  $  and  Flf,  10,  One 
recalls  that  in  tho  voltage  spectrum  there  are  throe  htoh  amplitude  reso¬ 
nant  peaks  at  f-1,  f50.8  and  f»0,375  and  a  number  of  smaller  peaks  at 
higher  frequencies,  In  Table  6  and  Fig,  10  one  observes  that* 

1  -  For  the  region  0,0< jx/d|<0,3,  only  the  resonance  at 
f«l  exists. 
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2  -  At  |x/dj  *>  0.4,  the  resonance  at  f»0.9  suddenly  appears 

and  the  resonance  at  £»C.575  is  noticeable.  The  resonance 

at  f*l  has  been  reduced  in  amplitude. 

3  -  At  jx/dj“0.5,  the  peak  f»l  has  been  greatly  reduced  while 

the  other  resonances  hove  been  reinforced. 

4  -  At  x/d-O.G,  the  resonance  at  £«1  does  not  show,  the  reso¬ 

nance  at  f*0.8  has  begun  to  decrease  in  magnitude  and  the 

lowest  resonance  has  become  predominant. 

5  -  At  x/d“0.7,  the  resonance  at  f«0.8  has  disappeared,  the 

resonance  at  f*l  shows  a  small  remnant, 

6  -  At  x/d*> 0.8,  only  the  resonance  at  approximately  f >*0.575 

remains. 

Also  worthy  of  note  is  the  high  magnitude  of  the  Fourier 
spectrum  at  x/c>0.3  for  frequencies  below  f»P.8,  This  tends  to  confirm 
the  "barrier"  behavior  discussed  nreviouslv.  It  io  also  possible  that 
some  small  resonances  are  hidden  in  this  part  of  tha  spectrum  which 
cannot  be  seen  due  to  the  limited  number  of  points  available. 

The  two  lower  frequencies  are  both  higher  than  the  local  cold 
plasma  frequencies. 

Thus  the  three  regions  are  established  and  can  be  defined: 

1  -  Region  £1,  in  the  canter  of  the  slab,  where  only  the  reso¬ 

nance  at  f*l  is  seen. 

2  -  region  #3,  near  the  walls  of  the  slab,  vrhere  the  resonance 

at  f-0.575  is  seen  and  only  traces  of  the  other  two  resonances 

are  seen  If  they  are  seen  at  all. 


'  "  n 

3  -  Region  i1 2,  between  the  previous  two  regions,  where  both  the 
previous  resonances  can  be  seen  plus  a  strong  resonance  at 
£«Q,8  which  is  characteristic  of  this  region. 

In  this  standard  experiment,  region  If  1  covers  0,0<jx/d|<0,3, 

Region  02  consists  of  the  two  sectors  0.4<jx/d|<0.6.  Region  if 3  is  formed 
by  the  two  sectors  defined  by  0.7<jx/di<0.8. 

The  fact  that  only  a  limited  number  of  points  per  frequency  unit 
are  available  might  of  course  hide  some  detail  fcut  this  will  not  influence 
the  above  results  to  a  significant  degree. 

C.  Effect  of  Variation  of  Parameters  of  the  Plasma  Model 

In  order  to  establish  that  the  two  lower  Frequency  peaks  are 
thermal  modus  as  conjectured  in  the  preceding  section,  the  effects  of  the 
variation  of  temperature  and  the  slab  width  are  investigated. 

Computer  runs  were  made  for  plasma  slabs  with  a  and  u  as  for 
the  standard  experiment  but  with  d«0.004 (0.0005)0.006  ,  0.007(0.001)0.010. 
Results  presented  here  are  for  d*=0.004,  0.005,  0.007  and  C.010  only  since 
these  suffice  to  describe  the  effect  of  variation  of  d. 

Figs,  11  and  12  show  the  voltage  as  a  function  of  tire  for  the 
slabs  of  halfwidth  d“0.Q04  and  d=0.010.  '..7:116  it  is  difficult  to  make  a 
comparison  between  Figs.  1  and  11,  one  is  able  to  notice  a  substantial 
difference  between  these  two  figures  and  ^ip.  1?..  In  the  latter,  the  peak- 
to-peak  times  are  much  more  uniform  than  in  the  former.  Peak  values  of  the 
voltage  and  times  of  occurrence  are  listed  for  all  four  widths  in  Table  3. 

"ore  interesting  and  inforrative  are  the  Fourier  transforms  of  the 
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voltages  for  the  four  plasma  3lahs  under  discussion.  Results  are  shown 
in  Table  5  and  plotted  in  Fig.  13.  One  notices  that 

1  -  The  peak  at  f^l  associated  with  the  central  region  (region  ill) 
increases  in  strength  with  increasing  slab  viidth,  but  with 
only  a  very  slight  shift  of  frequency,  indicating  a  cold 
oscillation. 

£  ■»  “he  next-peak  -associated- v; ith  the  "mid d  1  c ‘Ling I on  (region  II 3 ) 
decreases  in  strength  with  increasing  slab  width.  For 
da3.004,  0.^05  a  low  shift  is  clear  and  one  also  notices  a 
decrease  in  amplitude  of  the  peak  with  increasing  width.  For 
dB0.C07,  0.010  peaks  ate  lately  identifiable  and  no  clear  con¬ 
clusion  can  he  drti'.m. 

3  -  The  lowest  peak  decreases  in  strength  and  broadens  markedly 
as  the  slab  widens.  The  low  shifts  of  the  peak  towards  lower 
frequency  as  the  slab  widens,  are  clearly  seen. 

The  lot'  shift  of  the  two  lower  frequency  peaks  as  a  result  of  the 
widening  of  the  slab  indicates  the  thermo resonance  in  the  sense  of  Baldw  in2? 

op 

and  Ignat  ,  Or  the  thermotiodcs  associated  with  a  circular  Tonks-Dattner 
structure  are  the  acoustic  wave  resonances  bounded  between  the  sheath  and 
a  dense  pla'tta  region  beyond  which  the  low  frequency  oscillations  cannot 
penetrate. 

This  in  further  conf'rmed  bv  a  study  of  the  effect  of  varying  the 

background  temperature,  i.e.,  the  thermal  velocitv  a.  T'xperiments  were 
5 

performed  for  u«10  ,  d*0,0C5  -  the  same  stimulus  and  width  of  the  standard 


experiment  -  with  a-8-105,  9*105,  l.l'lO6,  1.2*106.  Results  of  the 
Fourier  transform  of  these  experiments  are  listed  in  Table  8,  together 
with  the  results  of  the  standard  experiment.  Fig.  15  plots  the  Fourier 
transform  for  the  highest  and  lowest  thermal  velocities  and  the  Fourier 
transform  of  the  standard  experiment.  One  notices*  by  comparison  with 
Table  5  and  fig.  13  that  a  decrease/increase  in  the  temperature  of  the 

equilibrium  slab  has  a  very  similar  effect  to  the  widening/narrowing _ 

of  the  plasma  slab  and  in  particular  that  the  highest  of  the  three  high 
amplitude  resonances  shows  a  negligible  frequency  shift  while  the  low 
shift  of  the  two  lower  peal'.s  is  again  clearly  noticeable. 

Thus  the  thermal  character  of  the  two  lower  high  amplitude 
resonances  as  well  as  the  cold  character  of  the  resonance  at  approxi¬ 
mately  f*l  has  been  clearly  established. 

The  Fourier  transform  of  AE(x/d,t)  for  d-O.GIO  was  also  obtained 
and  has  been  tabulated  in  Table  7  and  plotted  in  Fig.  14  for  x/d-0.0,0.5,0.8. 
The  behavior  is  in  accordance  with  the  result  on  the  voltage.  One  notes 
the  small  peak  around  f«0,7  for  the  transform  at  x/d*»0.8, 

D,  Simmarv  of  Results  in  the  Linear  Pcgion 

In  the  preceding  sections  it  has  been  shown  that  the  rlasma 
model  presented  has  a  linear  region  and  that  in  this  region  the  Fourier 
spectrum  of  the  impulse-response 

1  -  has  three  well-pronounced  resonance  peaks  below  and  around 
the  plasma  freoucncy  at  the  center  o^  the  slab,  and  that 
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FIG.  14  FOURIER  TRANSFORMS  OF  AE{t) 
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FIS.  15  FOURIER  TRANSFORMS  OF  V(t) 

a  =  8.105,  106,  1.2.1 06 ,  u  =  105j  d  *  0.005 
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2  •  each  of  these  freeuencies  is  associated  ”ith  a  relatively 
well-defined  re"ion  of  the  slab. 

It  has  als.,  been  she m  that  the  follovrinp  associations  can  -e 
made  between  Che  three  resonances  and -areas  of  the  slabs 

1  -  The  highest  of  the  three  resonances  belongs  to  a  region 

located  around  the  center  of  the  slab.  The  arrlitude  and 
frequency  arc  rclativcl”  independent  o£  plasma  temperature 
and  vidth  of  slal .  The  value  of  the  fre^ucnc”  of  the  resonance 
peak  is  at  or  slight!”  belo'-’  the  plasma  frequency  at  the 
center  of  the  slab, 

2  ~  The  lowest  of  the  three  resonances  is  associated  with  the 

region  near  ti.c  plasra  sheath.  The  frcnuenc”  of  this  reso¬ 
nance  peak  as  well  as  its  amplitude  as  ’.-ell  as  the  3-db  width  of 
the  resonance  depend  on  both  slab  width  and  temperature. 

The  frequency  of  the  peak  is  above  the  local  cold  plasma  fre¬ 
quencies  in  the  region, 

3  -  The  "middle"  resonance  is  associated  with  the  "middle"  region 

located  between  the  "center"  region  and  the  "wall"  region. 

The  amplitude  of  this  resonance  depends  on  slab  width  and 
temperature;  the  frequence  of  this  resonance  is  also  above, 
the  local  cold  plasma  frequencies  of  this  region. 

Thus  one  is  able  to  conclude  that 

1  -  The  resonance  peak  at  f*l  is  a  cold  plasma  resonance  associated 
with  the  plasma  frequency  in  the  center  of  the  plasma. 
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2  -  The  other  two  resonances  are  thermal  resonances. 

Table  9  shows  local  plasma  freauencles  and  periods  In  the  back¬ 
ground  plasma  for  x/d*0.0(0ll)0.9,  normalized  to  center  plasma  frequency 
and  period  of  unity. 

It  has  further  been  shown  that  the  three  regions  can  be  defined 
in  the  following  manner  through  the  Fourier  transform  of  AE(x/d,t) 

1  -  The  central  region  contains  only  the  resonance  at  f*l.  The 

other  resonances  cannot  be  seen  In  this  region  or  are  of  neg¬ 
ligible  magnitude. 

2  -  The  wall  region  shows  the  lowest  resonance  strongly  predomi¬ 

nant.  The  other  resonaitces  minht  or  might  not  be  noted  in 
this  region  and  are  of  small  magnitude  whenever  noticed. 

3  -  The  ’Viiddle"  region  contains  all  three  resonances  but  onl”  in 

this  region  is  the  "middle"  resonance  very  strong,  i.e.,  of 
the  order  or  larger  than  any  of  the  other  two  resonances. 
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TABLE  9 


Local  Plasma  Frequencies  and  Periods 
Normalized  to  1  at  Center  of  Slab 


X 

a)  ■  ■  or  f 

T 

d 

pn  -on 

pR 

0.0 

1, 

1. 

0.1 

.9877 

1.012 

0.2 

,9511 

1.051 

0,3 

.8919 

1.122 

0.4 

.8090 

1.236 

0.5 

.7071 

1.414 

0.6 

.5878 

1.701 

0.7 

.457:0 

2.203 

0.8 

.3030 

3.236 

0.9 

.1564 

6.394 
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V.  RESULTS  IS  THE  NON-LINEAR  REGION 

One  of  the  goals  of  the  present  investigation  was  to  device  a 
numerical  scheme  to  investigate  the  impulse-responses  of  a  non-uniform 
bounded  plasma  with  sheaths  under  strong  excitations  into  the  non-linear 
regime.  The  approach  presented  herein  fulfills  this  goal  to  some  extent. 
One -recalls  that  the  program  is  based  on  the  non-linear  Boltamann-Vlasov 
equation  with  the  imposed  Inhomogeneous  background  containing  the  singular 
fields  at  the  slab  boundaries  simulating  the  sheath.  The  limitation  im¬ 
posed  into  this  direction  comes  from  the  availability  of  computer  time  and 
of  computer  core.  This  has  already  been  discussed  in  Section  lll.b. 

The  time  during  which  the  solution  can  be  considered  valid  varies 
as  a  function  of  several  parameters.  One  important  factor  causing  the 
numerical  instability  is  the  singular .sheath  field.  Thus  in  a  different 
problem  such  as  that  of  a  unbounded  plasma  or  a  uniform  plasma  considered 
in  most  of  the  previous  investigations,  the  numerical  instability  problem 
should  be  considerably  relaxed.  In  fact,  it  is  observed  that,  if  the  slab 
is  widened,  the  situation  improves. 

With  the  computation  resources  available  to  the  present  study,  the 
strength  of  the  stimulus  has  been  extended  into  the  non-linear  regime  with 
some  limited  degree  of  success.  For  experiments  defined  by  a«10^  and 
d*>0.G05,  i.e.  the  temperature  and  slab  width  of  the  standard  experiment, 
the  validity  of  data  extends  to  about  7  plasma  periods  for  an  excitation 

e  c  g 

u*>2.10  ,  to  3  plasma  periods  for  u"»'5.10  and  to  2  plasma  periods  for  u-10 
These  times  can  be  extended  by  refining  the  grid  in  all  three  dimensions 
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but,  as  mentioned,  this  will  require  more  core  3pace  and  availability  of 

relatively  large  amounts  of  computer  time.  To  give  as  idea  os  the  volume 

£ 

of  calculations  involved,  in  the  standard  experiment  32.1CT  points  of  the 
electron  distribution  function  are  calculated  for  10  plasma  periods,  each 
one  requiring  several  additions,  multiplications  and  of  course  requiring 
several  data  which  have  to  be  brought  from  memory  co  the  arithmetic  proc¬ 
essor.  Further  developments  in  computer  technology  may  alleviate  the 
problem. 

From  the  computations  performed,  the  basic  wave  form  of  the 

voltage  remains  quite  similar  to  linear  results.  Some  of  the  obtained 

results  are  listed  in  Table  10  and  11.  This  listing  of  the  voltage  peaks 

and  their  time  of  occurrence  permits  to  inf*r  that  the  locations  of  the 

peaks,  at  least  during  the  relatively  early  stages,  will  not  change  to  a 

considerable  amount  but  that  the  amplitudes  at  the  peaks  might  change  to 

some  degree  which  at  the  moment  is  difficult  to  as certain. 

The  Fourier  transform  of  the  voltage  has  been  performed  and  results 

6  5 

are  given  in  Table  12  for  the  experiment  defined  by  a*»10  ,  u«*2.10  ,  d»0.005 
and  shown,  together  with  the  Fourier  transform  of  the  voltage  of  the 
standard  experiment,  in  Fig.  16.  Hie  qualitative  behavior. of  the  trans¬ 
forms  is  the  same.  A  difference  in  magnitudes  is  expected  due  to  the 
normalization  process  and  due  to  the  attenuation  expected;  one  Is  in  face 
of  the  same  phenomenon  discussed  for  results  of  Table  4  and  Fig.  3.  Briefly, 
one  notes  that: 

1.  The  resonance  of  the  peek  corresponding  to  the  region  near  the 

wall,  at  approximately  f«0.575,  shows  an  increase  of  approximately 
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TABLE  1£ 

\  '  ■ .  •  .  v 

.Envelop!  Points  of  the  Fourier  Transform' of  V(t) 
a  -  lO6,  u  *  2.1G5,  d  *  0.005, 


.000 

.637 

.013 

.674 

.265 

.7% 

.395 

.955 

.560 

2.2976 

t 

f 

.790 

1.7908 

.350* 

1.1364 

.985 

2.9311 

1.145 

.8480 

1.285 

.4974 

1.415 

.3430 

1.550 

.2759 

1.680 

.2154 

1.810 

.1866 

1.940 

.1494 

2.075 

.1382 

2.200 

.1128 

2.335 

.1118 

2.465 

.0859 

t  Results  indicate  probable  existence  of  phase  reversal 
point  between  0.650  and  0.590,  Exact  position  and 
absolute  value  at  that  point  cannot  be  determined. 

In  Figure  Lo.  16,  this  phase  reversal  point  is  ex¬ 
trapolated  with  a  certain  arbitrarity  to  show  the 
valley  between  the  two  resonances. 

*  Phase  reviatsal  point 


V(«) 
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2.  The  resonance  peak  at  approximately  0.8,  generated  in  the 
intermediate  region  (region  2)  shows  an  increase  by  about  6.5%. 

3.  The  resonance  corresponding  to  the  region  in  the  canter  of 
the  slab  is  approximately  31%  higher  chan  that  of  the  standard 
experiment. 

It  should  be  remembered  that  the  results  for  the  experiment 
under  discussion  had  to  be  obtained  from  a  time  domain  truncated  to  7.5 
plasma  periods.  Therefore  the  number  of  envelope  points  is  rather  small 
as  can  be  seen  from  Fig.  lo  and  Table  12.  The  adduced  results,  however, 
permit  to  affirm  that  one  has  begun  to  penetrate  the  non-linear  region. 

As  mentioned  before,  once  one  leaves  the  linear  region  An  does 
not  remain  identically  zero  at  the  center  of  the  slab.  Variations  of 
several  percent  of  equilibrium  density  at  the  center  of  the  slab  are 
observed  for  a“u«*10^,  d»0.CQ5.  The  results  are  of  the  same  magnitude 
if  d-0.010.  This  implies  that  the  particle  migration  pattern  ar  the 
origin  will  no  longer  have  the  "dipole  characteristic"  around  the  center 
of  the  slab.  Further,  Aa  will  no  longer  be  essentially  odd  nor  will  AE 
be  essentially  even. 

The  non-linear  interaction  for  the  high-amplitude  excitation  is 
vividly  shown  in  this  study:  Consider  the  case  a^u^lO^,  d=0.005,  i.e. 
the  configuration  or  the  standard  experiment  with  a  ten  times  stronger 
stimulus.  It  is  possible  to  see  that  a  strong  tendency  to  develop  bunch 
ing  is  taking  place  at  the  very  early  stages.  In  Fig.  17,  the  total 
electric  field  (not  its  deviation)  is  plotted  as  a  function  of  x/d.  One 
observes  initially  the  field  at  t-0.3  (at  about  the  first  voltage  peak) 


93 

which  is  characteristic  for  the  electric  field  during  the  second  quarter 
of  the  first  placaia  period,  i.e.  0.25<t<0.5.  The  electric  field  is 
negative  everywhere  except  very  near  the  wall  x«d;  practically  all 
electrons  which  were  accelerated  by  the  initial  impulse  towards  the 
wall  at  x»-d,  now  she  a  force  which  tries  to  restore  them  towards  the 
wall  at  x«d,  This  is  similar  to  the  linear  case.  Next  one  observes  the 
field  at  f*0.6  (at  about  the  first  voltage  zero)  which  is  characteristic 
for  the  electric  field  duritjg  the  third  quarter  of  the  first  plasma 
period,  to  be  more  exact,  fer  thfe  time  0.55<t<0.70.  One  sees  a  positive 
region  in  the  middle  and  a  r.qgative  region  between  the  center  of  the  slab 
and  the  wall  at  x“d.  One  observes  three  zeros  of  the  electric  field  while 
in- the  linear  region  one  observes  only  one.  The  particle-wave  interaction 
is  not;  seen  clearly;  the  electrons  In  the  regions  of  positive  electric 
field  are  accelerated  towards  the  "left",  i.e.  the  wall  at  x»-d  while 
electrons  in  "the  regions  of  negative  electric  field  are  accelerated  towards 
the  "right",  i.e.,  towards  the  wall  st  x»+d.  Thus  bunching,  is  produced 
around  x/d*-0,35  and  some  around  z/dS+0.75  while  a  depletion  takes  place 
around  x/u*+o.ty_>. 

Also  shown  in  Fig.  17  fer  comparison  purposes  is  the  electric 
field  in  the  plasma  in  equilibrium. 

Similar  buncoing  and  depletion  areas  are  also  observed  at  later 
times  during  the  experiment.  Similar  results  are  observed  for  a«u“10b, 


d-0.010. 
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APPENDIX  Ii  TliL  GENERAL  SYSTEM  GF  EQUATIONS  FOR  A  COLLISIONLESS  GASEOUS 

PLASMA.  REDUCTION  TO  THE  MONODIMENSIOKAL  (LONGITUDINAL)  CASE. 


The  system  of  equations  for  a  collisionless  gaseous  plasma  formed  by 
n  kinds  of  particles  is 


3E 

tvr  ~  V  x  11  +  /  f  vdv 

oct  ~  a  ot  *■ 


(l.l) 


3 II 

P  x  g  +  yo~ 


a. 2) 


3f  q 

~+v?f  +  — (E+y  vxH)  •  V  £  «  0 
3t  a  nu,  ~  o~  ~  v  a 


(1.3) 


7  •  5  "  Zaf  !  V* 

o 


(I. A) 


v  .  h  -  o<-  y.u) 


(1.5) 


B  -  y  5i 
~  o~ 


(I.C) 


U  *  1  ,  2  ,  3  ,  » •  • ,  n 


The  8+n  scalar  equations  (1.1-5)  relate  the  6+n  scalar  functions  E,H,fa  to 
the  7  independent  variables  £,v,t.  ^*e  tYPe  particle  is  characterized 

by  the  subscript  a.  e0  and  yQ  are  the  permeability  and  permittivity  of 
free  space.  E  and  h  are  the  usual  field  vectors,  B  being  defined  by  (1.6) 
and  the  respective  result  used  in  (1.3).  f  is  the  distribution  function 
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of  the  a-th  type  of  particle.  The  system  (1.1-6)  refers  the  variables  within 
the  international  M.K.S.  (Giorgi)  or  similarly  rationalized  system. 

Initially  assume  that  one  has  a  gaseous  plasma  composed  of  one  ele¬ 
ment,  the  amount  of  impurities  being  negligible.  This  gas  is  fully  ionized. 
It  Is  not  necessary  to  assume  that  the  plasma  is  compensated.  One  has  only 
electrons  and  ions;  n*2.  Assume  that  one  is  interested  in  a  relatively  short 
time  span  and  that  initially  the  ion  velocity  can  ba  considered  small;  then 
the  ions  can  be  considered  immobile.  Thus  the  ion  distribution  function  does 
not  depend  on  velocity  or  time  but  on  space  only. 

Consider  now  a  plasma  contained  by  the  two  walls  x*  +<3.  If  one  im¬ 
poses  translational  invariance  with  respect  to  y  and  z  and  further  imposes 
that  the  electrons  move  in  the  x-direction  only,  the  electron  distribution 
function  will  be  function  of  the  space  coordinate  x,  the  velocity  v  (  in  the 
x-direction)  and  time  t  only.  From  (1.4)  follows  that  aiv  E  is  a  function 

of  x  and  t  only;  assume  the  same  to  be  true  of  the  component  E  which  one 

X 

sets  £*2^  for  notational  simplicity.  One  notes  further  that  (1.3)  reduces  to 
0*0  for  the  ions  and  that  the  dot-cross  product  of  the  same  equation  for 
electrons  vanished  since  the  velocity  gradient  is  collinear  with  v  both  being 
in  the  x-direction.  Thus  one  lias 


^  e 

dX  "  ~  £ 

I  f  dv  +  |—  n 

(1.8) 

0 

0 

3f  3f 

?t  'cX 

e  a 

—  2x—  “  0 
m  9v 

(1.9) 

f  being  the  electron  distribution  function,  n(x)  the  ion  density  and  e  now  is 
the  absolute  value  of  the  electron  charge  and  It  is  assumed  that  the  charge  of 
an  ion  is  e. 


It  should  bs  mentioned  onae  mere  Chat  the  plasma  does  not  have  to  be 
compensated;  in  fact  it  ie  more  convenient  in  further  development  to  consider 
the  plasma  not  compensated  since  this  will  lead  to  ease  of  language.  One 
notes  that  (1=9)  contains  no  inform cion  about  the  ion#  except  through  £ • 

Xaus  the  electron  distribution  function  will  not  be  affected  if  the  term  in 
n(x)  be  replaced  by  the  x-derivative  of  a  time-invariant  external  field.  It 
is  of  course  also  possible  that  this  external  electric  field  has  the  x- 
derivative  such  that  the  plasma  becomes  compensated. 

With  these  assumptions  one  recognizes  in  the  system  (1.8,9)  equations 
for  a  longitudinal  plasma  behavior.  It  will  be  shown  that  (1.8,9)  are  in 
fact  sufficient  together  with  mixed  initial  value  and  boundary  conditions 


to  determine  f  and  £. 


91 


APPENDIX  It:  PP.VT  OF  "  E  SEFAPA3ILITY  OF  VAR1ABTE8  OF  THE  '  IST'irUTTOU 
FUNCTION  f(x,v)  OF  THE  CONSIDERED  PLASMA  SLAB  If!  ITS 

eqtiii.ibp.tuii  state. 


In  this  appendix  we  prove  that  for  E(x)  real  analytic,  the  distribu¬ 
tion  function  f(x,v)  is  separable,  i.e.,  that  f(x,v)  "^(x)iJj(v)  and  that, 
further,  ii>(v)  is  Gaussian  and  that  either ? (x)  >  0  ory(x)  <  0  everywhere 
in  the  slab. 

Ne  consider  the  following  equilibrium  state  of  a  Vlasov  plasma 
+oo 

-  -  /  f(x,v)dv  +  n(x)  (II. 1) 

-00 


.  EU) 


defined  on  the  domain  D 

-  d  <  x  <  d 


oo  <  v  <  oo 


and  subject  to  the  boundary  conditions 

f(-d,v)  «  f(d,v)  »  f(x,-co)  «=  f (x ,+co)  ■  0 
C  is  a  real  positive  constant. 


(II. 2) 


(II. 3) 


(II. 4) 


Lemma:  diven  f(x,v)  differentiable  everywhere  or.  D,  E(x)  and  n(x)  are 
uniouely  determined  provided  dE/dx  exists  for  all  xeD  and  the 
integral  in  II. 1  exists. 

Proof:  f(x,v)  is  differentiable  everywhere,  the  integral  in 

(11,1)  exists.  The  rest  follows  by  inspection  from  (11.1,2), 


1 
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Lamina:  If  f(x,v)  is  differentiable  and  separable  in  its  variables,  then, 
k  and  a  being  constants 


f(x,v)  •  fc<p(x)exp(-vz/2az) 


(11,5) 


£  : 


Proof :  Let  f(x,v)  *  <P(x)iHv)  •  Then  (71,3)  yields 


\ 


E 


£*£ 

:js  v& 

Cj)  d4> 

dv 


(11,6) 


Since  E  *  E(x),  the  second  fraction  has  to  be  a  constant.  Integrating  the 
resulting  ordinary  differential  equation  v?ith  K  and  a  constant 


£(v)  *  K  exp(-v*/2az) 

and  (II. 5)  results.  It  is  convenient  to  choose  K  such  that 

-co 

/  v(v)  dv  •  1 

—CO 

resulting  in  the  Gaussian 

9(v)  m  — 1—  exp(-vx/2a2) 
av'H 


It  is  also  convenient  to  cheese 


o(0)  *  1 


(II. 7) 


(II. 8) 


(11.9) 


(11.10) 


Lemma:  If  f(x,v)  is  separ-ble,  i.e.,  f(x,v)  »tp(jt)'Mv),  then 

F(*)  -  -  •  -  a2v— (fn<p) 


a2  d2 

lux)  »  y  -  ^ - -(Lnc.) 

dx 


provided  the  derivatives  in  (11.11,12)  exist. 


(11. 11) 

(11.12) 
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Proof;  One  obtains  E(x)  be  er.teritig,  (II. 5)  into  (11,2).  Then,  by  entering 
(11.5,11)  Into  <11. 1),  (11.12) fellows. 

THEOREM !  Let  £<x)  be  a  real  analytic  function  in  0  and  let  v(x)  be  a  solution 
of  the  ordinary  differential  equation  (11.11) 

jjjttop)  »  -  Kny.)  (11.13) 

satisfying  9(-<i)  "  o(d)  ®  0.  Then  f(x»v)  *  K9(x)  axp(-v2/2az) , 
i.e.,  f(x,v)  is  separable  in  its  variables  and  its  velocity 
dependence  is  Pauseiuo. 

Proof  >  Let  f(x,v)  =  '+(:<} \p(v)  where  iHv)  is  given  Jy  (11.9)  and 

V  (x)  is  the  solution  of  (11.13)  satisfying  <p  (d;*cp(-d)  »  0.  f(x,v) 
together  with  the  viven  E(x)  satisfies  (II, 2)  as  well  as  the  boundary 
conditions  (II. A).  Eut  for  'E(x)  analytic,  the  theorem  of  existence 
end  nnicityof  the  solution  of  cuasi-linear  partial  differential  equations 
of  the  first  order  in  two  variables  certainly  applies.  Thus  F(x,v)  * 
<?(x)’y(v)  will  also  be  the  only  solution  of  (II. 2)  with  the  given 
E(x) .  O.E.D. 

THEOREM:  If  E(x)  is  analytic  in  P,  9  is  either  positive  or  negative  everywhere 
in  P. 

Proof :  Prom  (iLll)  follows,  using  (II. 10),  that 

9(x)  =*  K  exp(-  j  /*  E(x)dx  )  (11.14) 


ICO 


Since  E(x)  is  analytic,  it  is  continuous  and  thus  <?(x)  is  continuous.  But 
E(x)  also  is  finite  lr,  D  anu  thus  tp(x)  is  con-aero  in  ?'  as  follows  from 
(XT. 11).  Therefore  <P(x)  is  either  always  positive  or  negative.  Q.E.D. 

One  notes  that  if  (11,3.0)  is  assured,  cp(x)  is  positive  every¬ 
where  in  D. 
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